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ABSTRACT 

The Carnegie Hubble Program (CHP) is designed to calibrate the extragalac- 
tic distance scale using data from the post- cryogenic era of the Spitzer Space 
Telescope. The ultimate goal of the CHP is a systematic improvement in the 
distance scale leading to a determination of the Hubble Constant to within an 
accuracy of 2%. This paper focuses on the measurement and calibration of the 
Galactic Cepheid Period-Luminosity (Leavitt) Relation using the warm Spitzer 
IRAC 1 and 2 bands at 3.6 and 4.5 /mi. We present photometric measurements 
covering the period range 4-70 days for 37 Galactic Cepheids. Data at 24 phase 
points were collected for each star. 

Three PL relations of the form M = a(log(P) — 1) + b are derived. The 
method adopted here takes the slope a to be -3.31, as determined from the 
Spitzer LMC data of Scowcroft et al. (2012). Using the geometric HST guide- 
star distances to ten Galactic Cepheids we find a calibrated 3.6 /tm PL zero-point 
of —5.80 ± 0.03. Together with our value for the LMC zero-point we determine 
a reddening-corrected distance modulus of 18.48 ± 0.04 mag to the LMC. 

The mid-IR Period-Color diagram and the [3.6] — [4.5] color variation with 
phase are interpreted in terms of CO absorption at 4.5 /mi. This situation com- 
promises the use of the 4.5 /im data for distance determinations. 

Subject headings: Cepheids distance scale infrared: stars Galaxy 
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Introduction 



The Carnegie Hubble Program (CHP) is designed to reduce systematic uncertainties 
in the distance s cale. The compelling reasons for doing so are provided in an overview by 
Freedman et al.l (120 111 ). The first phase of the CHP is a warm Spitzer legacy mission, the 
preliminary goal of which is to reduce the systematic uncertainties in H to 3% or better. 
The second phase will include observations from GAIA and the James Webb Space Telescope 
(JWST) where the goal will be to push this number to 2%. 

The warm Spitzer phase consists of observing Cepheids in the Milky Way (MW), Large 
Ma gellanic Cloud (LMC), Small Magellanic Cloud (SMC) and other Local Group Galaxies 



see 



Freedman et al.ll201ll for a complete list) to calibrate the local distance scale in the mid- 



infrared (mid-IR). The program extends into the Hubble Flow by calibrating the mid-IR 
Tully-Fisher relation and fart her to Type la supernova host galaxi es observed as part of the 
Carnegie Supernova Project ( IFolatellill2009t IContreras et al.ll2010f ). 



This paper presents warm Spitzer IRAC channel 1 (3.6 /im) and channel 2 (4.5 /im) 
light curves for 37 Galactic Cepheids. Each star was observed at 24 phase points. The 
data are used to derive robust mean magnitudes and colors, and, with the adoption of HST 
parallaxes for 10 stars, an accurate calibration in the mid-IR. The absolute magnitudes 
give a distance to the LMC, which we believe is currently the value carrying the lowest 
systematic uncertainty. The calibrations and LMC distance values for a sample of Cepheids 
in Galactic Clusters and for a sample with Infrared Surface Brightness data are compared 
to the HST-based distance. A brief discussion of the observed Period-Color relationship and 
its interpretation in terms of CO affecting the 4.5 /im band is also given. 



2. Warm Spitzer Observations 
2.1. Target Selection 



The 37 Cepheids in our sample h ave multiple distance and extinction estimates, and span 



a wide range of 4 - 7 days in period. ( IFernie et al.lll995l . iTammann. Sandage. fc Reindlll2003 



Fouque et al.ll2007l ). The sample includes the majority of the nearest Cepheids; these should 

soon have high-precision parallaxes available from the GAIA mission (jWindmark. Lindegren. &: Hobbs 
20111 ). Table [TJ lists the target star, the adopted reddening and three measurements of dis- 
tance for each star if available: direct parallax data from HST, distance moduli obtained 
from main-sequence fitting of their host clusters, and distances from Infrared Surface Bright- 
ness (IRSB) measurements; the latter two methods were converted to units of parallax for 
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comparison. Reddening estimates from photometric and spectroscopic methods are listed in 
Table [2] as well as space reddenings determined from stars along the same line of sight. The 
average reddening of these methods was used in Table [TJ 



2.2. Observations 

Observations were made using the Spitzer Space Telescop e as part of a two-ye ar Ex- 



ploration Science Program, PID 60010: The Hubble Constant (IFreedman et al.ll2008l ). The 
warm Spitzer Mission started in 2009 (Cycle-6) and the Galactic Cepheid observations were 
completed in early 2011. Each Cepheid was observed at 24 epochs, pre-selected and sched- 
uled to fully sample the light curve (23 epochs for ( Gem). At each epoch a nine-point dither 
pattern was used to mitigate array- dependent artifacts such as bad pixels and cosmic rays. 

The majority of the data were taken using the sub-array mode of the Infrared Array 
Camera (IRAC0) ( Fazio et al.l 120041 ). with the shortest available frame time of 0.02 seconds 



(effective exposure time of 0.01 s). The sub-array mode outputs data from only one corner 
of the detector, in a 32 x 32 pixel format, thus allowing the shortest possible exposure times. 
The observations of CF Cas were made using the full array mode (0.4 s frame time, 0.2 s 
effective exposure time) because it is relatively faint compared to the other program stars. 

The sub-array data are provided by the Spitzer Science Center (SSC) in two forms: as 
an image cube of 64 frames (each 32 x 32 pixels) and as a single combined image; all further 
discussions to sub-array data refer to the combined form (sub2d extensions). All data were 
retrieved in the basic calibrated data (BCD) format, and were reduced using the most recent 
pipeline (S18.18.0). 



3. Data Reduction and Photometry 



The stellar flux in each i mage was measured using the Mosaicking and Point-source Ex- 
traction (MOPEX@) package (IMakovoz fc Marleaull2005l ) . The pipeline script apex_userJist-lframe.pl 
(included with MOPEX ) was used to perform profile fitting photometry. Uncertainty images 



lr The IRAC instrument handbook and ancillary data products are available at: 
|http://irsa.ipac.caltech.edu /data/SPITZER/ docs/irac/| 

2 The MOPEX software and documentation are available at: 
|http://irsa.ipac. caltech.edu /data/SPITZER/docs/dataanalysistools/tools/mopex/ 
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were supplied to MOPEX using the square root of the input image fram^. 



Once the first-pass fluxes had been obtained we processed the data to correct for three 
systematic effects. These are (1) the masking of saturated or markedly non-linear pixels; 
(2) corrections for non-uniformity of response across individual pixels; and (3) correction 
for image persistence, in which a bright source will leave behind a trail of spurious flux as 
the telescope executes its dither pattern. The three effects are examined and the derived 
corrections are explained in detail in Appendix |A] 

The final stellar flux for each epoch is determined from the mean of the 9 dithered 
flux measurements from MOPEX, modified by the three corrections. The random error is 
adopted as the dispersion of the 9 measurements and the systematic error is taken as the 
zero- point error adopted by the SSC, viz., 0.016 mag for both [3.6] and [4.5] ( jReach et al. 
20051 ). Table [3] shows a sample of the IRAC photometric data (the magnitudes are named 
[3.6] and [4.5]) available for the 37 starfl 



Results 



4.1. Mid-IR Light Curves for 37 Galactic Cepheids 



Periods from the General Catalog of Variable Stars (GCVS, ISamus et al.l 120091 ) were 
assumed. In the cases of U Car and V340 Nor, periods were computed using photometric 
data from lLaney fc Stobid (Il992f ) . Figure [T] presents the individ ual light and c o lor cu rves 
(Vega magnitudes) for each Galactic Cepheid. Data points from iMarengo et al.l (120101 ) are 
shown as open triangles for comparison, when available. All are in good agreement except 
for U Car w here the difference is lik ely due to a phase shift resulting from a period increase 
between the lLaney fc Stobid (119921 ) data and ours. All the light curves are plotted with the 
same magnitude range to emphasize relative changes in signal-to-noise ratio and amplitudes. 
The internal photometric precision is high, ranging from 0.004 to 0.029 mag. Interestingly, 
for the longer-period Cepheids, there is strong variation in the [3.6] — [4.5] color. This is due 
to the temperature-dependent carbon monoxide (CO) band in the [4.5] Channel, as will be 
discussed in 

Smooth light curves were generated using a Gaussian local estimation (GLOESS) al- 



3 The input values were converted to units of electrons first and the uncertainty images were converted 
back to units of MJy sr _1 . 

4 The full table can be accessed in the online version of this paper. 
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gorithm. GLOESS is an interpolating method that uses second-order polynomials to fit 
the data locally throughout the cycle. The data points surrounding the point to be fit are 
assigned weights according to a Gaussian window function; weights depend on their dis- 
tance (in phase ) from the fit point. The method has been used to conveniently obtain mean 
magnitudes by iPersson et al.l ( 12004 ) and IScowcroft et al.l ( 120111 ) for LMC Cepheids, and by 



Monson fc Pierce! (120111 ) for Galactic Cepheids. These data are uniformly sampled so the 
error on the mean is: a = tt4t^, where A is the peak-to-peak amplitude o f the light curve 



JV\/i2' 



and N is the number of sample points; this is discussed in the Appendix of IScowcroft et al. 
(120 111 ). The final total uncertainty in the mean magnitude turns out to be dominated by 
the systematic zero-point calibration of the Spitzer warm mission. Because the uncertain- 
ties in Channels 1 and 2 are correlated with each other, uncertainties in the color were also 
determined using the above equation. Table H] gives the [3.6] and [4.5] IRAC intensity-mean 
magnitudes and colors for the 37 Cepheids. 



4.2. Mid-IR Extinction Corrections 



Before discussing the Period-Luminosity and Period-Color Relationships, we shall need 
to correct for extinction. Compared to optical wavelengths, reddening and extinction cor- 
rections are relatively small at mid-IR wavelengths. They must, nevertheless, be quantified 
and applied because values of Ay can exceed 3 mag in our sample. We adopted an ex- 
tinction law for all stars that is applicable a long average lines - of-sigh t through the diffuse 
in terstellar medi u m. T he extinction law of llndebetouw et al.l ( 120051 ) combined with that 
of ICardelli et al.l ( 119891 ) best fulfill this choice. The relations: A\^^/ A K = 0.56 ± 0.06 and 
A[4.s\/Ak = 0.43 ± 0.08 (llndebetouw et al.ll2005l ) were derived from field stars in the Galactic 
Plane and are probably applicable to the Cepheids in this stud)|E 

To scale the extinctions at K to the reddenings E(B — V) we used the extinction law de- 
rived by ICardellLelaD ([ljfl): A(X)/A V = a(x) + b(x)/R v , where a = 0.574a: 1 - 61 , b = -0.527a; 1 - 61 , 
x — 1/A and Ry is the ratio of total-to-selective absorption (R v = A V /E( B — V)). Using a 
wavele ngth of A = 2.164 /zm for the K filter (actually K s ) as adopted by llndebetouw et al. 
(120051 ) and an average Ry = 3.1, we find: Ax/Ay = 0.117. The combined relations yield a 
final total-to-selective extinction of A [3 _ e] /E[B - V] = 0.203, A [4 . 5] /E[B - V] = 0.156, and 
J5([3.6] - [A.h])/E[B — V] = 0.047. 



Use of the llndebetouw et al.l ( 120051 ) mid-IR extinction law might be questioned on gen- 
eral grounds. Three stars have E(B — V) > 1, above which the corrections will begin 



5 Following 



Indebetouw et al. 



(|2005l ) we use K to mean the K s filter of the 2MASS survey. 
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to i ntroduce systematic e rrors. Fo r example, measured values of A^a/Ax range from 
0.4ldChapman et all l2009h to 0.64 rtFlahertv et aD 12007^ . Toward the Galactic Center 
Nishivama et all ( 12009^ obtain 0.50 ± 0.01. For V367 Set (E(B - V) = 1.231) the total 



range in v4[ 3 6 ] is 0.10 mag. The uncertainties in E(B — V) are all < 0.03 mag except for 
for GY Sge, where a(E(B — V)) = 0.17. This value introduces uncertainties of 0.035 and 
0.027 mag into the 3.6 and 4.5 fim magnitudes, respectively. Finally, if the uncertainty in 
A[ 3 ,q]/E[B — V] is as large as 0.05 (see above), the corresponding corrections will remain 
negligibly small. 



5. Period-Luminosity Relations at 3.6 and 4.5 /zm 

We now present the Period-Luminosity Relations for the 37 Cepheids in our sample. 
Table [1] shows that the sample may be divided into three subsamples, depending on the 
origin of their distance measurements. Henceforth we consider each subsample separately, as 
the three methodologies for distance determinations are quite different. In §5.21 we discuss 
in detail three weighting techniques for the data points. In the following sections we have 
adopted unweighted fits in finding slopes and zero-points. 



5.1. The Three Subsamples 



Ten of the 37 Cepheids have dire ct geometric parallax es determined from Hubble Space 
Telescope (HST) guide camera data ([Benedict et al.l 120071 ) and therefore have the most ac- 
curate distance determinations currently available. The data for the sample are provided in 
Table [51 Figure [2] shows the data and zero-point fit using uniform weighting; the distance 
uncertainties are displayed for reference as error bars. The data include the final de-reddened 



[3.6] and [4.5] magn itudes, the fina 



nitudes. Following 


Benedict et al. 


(Lutz & Kelker 


1973; 


Hanson 


1979 



adopt ed distance moduli, extinctions, and absolute mag- 



(120071 ) we have applied Lutz-Kelker-Hanson corrections 
LKH) to these parallaxes; the corrections are systematic 
and range from -0.02 to -0.15 mag with uncertainties of ±0.01 mag. For completeness we 
include in Table [6] and Figure [3] the corresponding data for the HST sample without LKH 



corrections. 



Eighteen of the 37 Cepheids are likely to be members of star clusters or associa- 
tions for which distance moduli have been estimated from Main-Sequence (MS) fitting 



6 See www.pas.rochester.edu/ emamajek/memo _ir_reddening.html for a summary. 
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dTurner fc BuAel 12OO2I : iTurnerl bold : iMaiaess et all l2012bl lal boilh . Table [7] lists the data 
and Figure H] shows the forced LMC-slope and zero-point fit for the uniformly weighted data. 
The Cepheids CEa and CEb Cas are presumably at the same distance as CF Cas by virtue 
of common membership in the cluster NGC 7790 and although separated by only l'.'O they 
were easily split using PRF photometry and they have been included in the sample. 

T hirty-two of the 37 Cepheids have distance determinations based on the IRSB tech- 
nique (IStorm et al.ll2011bl . and references therein). Table |H] contains the data and Figure 
shows the fit and zero-point for uniformly weighted data; W Sgr was rejected from the 
analysis due to its relatively high uncertainty. 

Each of the the data subsamples were fit using a PL relationship of the form: M = 
a (log P — 1) + b. As will be shown below we find no statistically significant difference between 
the slope of the [3.6] PL of —3.31 ±0.05 for the LMC, and that of the HST parallax sample. 
We thus adopt the LMC slope for the PL fit and re-determine zero-points for each subsample. 
The magnitude residuals from the PL fits are highly correlated with each other, suggesting 
that the widths of the PL relations are not driven by random photometric errors. Rather, the 
correlated scatter is most likely some combination of deterministically correlated (unit slope) 
distance errors and the intrinsic (correlated) positions of these Cepheids in the instability 
strip (IS). If the IS is represented by a rectangular distribution, (i.e., it is uniformly filled and 
has hard limits at the blue and red edges) then the peak-to-peak width in the residuals can 
be interpreted as the width of the IS or at least an upper limit, which in the HST subsample 
is ~ 0.4 mag. 



5.2. Dependence of PL Relations and Uncertainties on Weighting Techniques 

The PL relations were fit to each of the three subsamples using multiple weighting 
schemes and also by further restricting the subsamples by period cuts. 

The final uncertainty in absolute magnitude for an individual Cepheid is dominated by 
the uncertainty in its distance. In deriving a PL relation, however, an additional spread is 
caused by the finite width of the IS, and biases in the PL slope and zero-point may result 
depending on how the strip is filled. To investigate these uncertainties and their effects on 
the derived PL fits, we applied different weighting schemes. In addition, for each of the data 
sets, we investigated different period cuts so that the Galaxy data sets more closely matched 
the period range of the LMC sample, viz., 6-60 days. Finally, a fixed slope determined 
from the LMC data was force fitted to the Galactic data to determine only the Galactic 
zero-point. The data were fit using a PL relationship of the form: M = a(logP — 1) + b. 
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Table |9] presents the results using three different weighting methods for the data in this 
analysis. 

The first weighting method applies a uniform uncertainty of 0.1 mag to each Cepheid, 
the purpose of which is to provide results presumably less biased by Cepheids which may 
have underestimated distance uncertainties. The second method applies traditional weights 
as o~~ 2 to the absolute magnitudes. The third method falls between the first two in that it 
assumes an intrinsic scatter in the IS. In this case an additional uncertainty of 0.1 mag is 
added in quadrature to each individual uncertainty. The value of 0.1 mag is adopted from 
the average RMS scatter of the LMC data points around the best fit. 



5.3. PL-Slopes 

The IRSB slopes closely match the HST slopes, which is not surprising because the most 
recent IRSB distances used a projection factor (p-factor) calibrated using the HST parallaxes 
(but without LKH corrections) as priors. The IRSB slope is better constrained because of 
sample size, but is still dependent on the adopted p-factor. The effect of the varying the 
weighting method is most noticeable in the Cluster Main-Sequence (MS) fits where the slopes 
differ by more than 2a between the first two methods, and converges to within l a of the 



HST, IRSB and LMC slopes using the third weighting method. As discussed by iTurner 



(120101 ) some long-period Cepheids populate the blue edge of the IS and can bias the slope, 
so it is necessary to include an estimate of the intrinsic width of the IS to reduce the bias. 
The benefit of the third weighting method is that the intrinsic width of the IS is included as 
well as individual uncertainties for each Cepheid. As can be seen in Table [HI the slopes for all 
three methods agree very well with each other using the softened weights. Since the slope 
from each method agrees with that of the LMC, we chose to adopt the better-determined 
LMC slope and to redetermine the zero-points for both the [3.6] and [4.5] PL relations. This 
decision is further backed by r ecent studies that fi nd near-identical PL slopes for the MW 



and LMC in the near-infrared (jStorm et al.ll2011al ). 



5.4. PL-Intercepts 

As mentioned above, some of the long-period Cepheids occupy the blue edge of the IS 
and although we forced a fixed (LMC) slope, the zero-point can now be slightly biased if 
the entire sample does not uniformly populate the IS. We therefore limit ourselves to adopt 
zero-points from the uniformly weighted fits. This effectively assumes the width of the IS is 
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the only source of uncertainty and can be treated as equal for each Cepheid. We also chose to 
make use of the entire period range which provides a larger sample and will more uniformly 
populate the IS. With these choices in hand, we now have zero-points of: HST (with LKH) 
= -5.80 ± 0.03, HST (without LKH) = -5.74 ± 0.03 , MS = -5.75 ± 0.05 and IRSB = 
—5.74 ± 0.02. We notice again the agreement between the HST (without LKH) and IRSB 
zero-points as must be the case (see above). The average LKH correction is -0.06 mag, which 
if applied to the IRSB calibration would shift the [3.6] IRSB PL zero-point to —5.80 ± 0.03 
mag. Because they are calibrated using the HST parallaxes the IRSB zero-point does not 
offer an independent baseline measurement, however it does better sample the instability 
strip and since the it yields the same zero-point and scatter it indicates that the HST data 
is not too effected by paucity. 

The Cluster Cepheids do offer an independent check of the zero-point and they appear 
to confirm the HST (without LKH) zero-point. We note, however, that the outliers S Vul 
and TW Nor are more than 10% discrepant compared to their IRSB distances and rejecting 
them would change the Cluster z ero-point to -5.79 mag; in agreement with LKH. This is in 



agreement with the iNgeowl ( 120121 ) results wh ich also find a 0.04 m ag relative offset b e tween 



their derived Wesenheit PL distances to the Storm et al. ( 2011b ) IRSB and Turner (hoioh 



1973. 


Hanson 


1979, 


Smith 


2003, 


Loredo 


2007) 



the LKH correction and an uncertainty in this correction of 0.01 mag ([Benedict et al.ll2007l ). 



We note that the scatter in the HST data is less than the average uncertainty assigned 
to the HST parallaxes and that we have chosen to adopt this (smaller) empirical scatter as 
a measure of the total zero-point uncertainty. 



The Distance to the LMC 



As will be discussed in §7] the 4.5 /j m data are likely to be affected by CO absorption 
while the 3.6 /im data are not (see also iFreedman et al.l ( 120121 )). Consequently, we have 
adopted the absolute PL zero-point from the HST Leavitt Law at 3.6 /im (—5.80 ±0.03 mag) 
and compare that with the apparent zero-point of the LMC PL relationship (12.70 ± 0.02 
mag), both zero-points are defined at logP = 1.0 and the PL relations are parallel to each 
other. We found no measurable metallicity effects for the MW and LMC C epheids at 3.6 
microns for which there are spectroscopic [Fe/H] values (IFreedman et al.ll2012l. Fig ure 2.). By 



adopting a net extinction to the LMC of E(B-V)=0.1 (IFreedman fc Madoreil2010l ) and using 
the extinction law discussed in §4.21 (which yields a total LMC extinction of Ap.gi = 0.02 
mag), we find a distance modulus (m — M)[ 36 i for the LMC of 18.48 ± 0.04 mag. 
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Alternatively, we followed a multi-wavelength approac h to solve for reddening and dis- 
tance modulus simultaneously (IFreedman &: Madord 120101) . The multi- wavelength photom- 
etry for the LMC are taken from: lUdalski et al.lll999l (B, V and 1^ ), iPersson et al.l 12004 



(J, H and K s ) and IScowcroft et al 



were compiled from the literature ( 



20121 f[3.6] and [4.5]) . The Galactic Cepheid photometry 



BerdnikoyJl2008l Barnes et al.lll997l : lLanev & Stobidll992 



Welch et al.lll984t iMonson fc Piercdl201ll ) and average magnitudes were found in the same 



manner as discussed in I4.1( see Table ITOl The multi- wavelength PL relations are shown in Fig 
Hand summarized in Table [TT1 The slopes were found by fitting the LMC data in the period 
range 3.8-60 days; a universal slope is considered. The Galactic zero-points were found using 
the LMC slopes at each wavelength using only the HST parallax Cepheids; see 15.41 for de- 
tails. The apparent distance moduli are plotted ag ainst inverse wavelength in Figure |HJ The 
standard extinction law from ICardelli et al.l (119891 ) was fit to the data to find the true LMC 
distance modulus of 18.48 ± 0.03 and average LMC reddening of E(B - V) = 0.12 ± 0.01 
mag. The K s and [4.5] data were excluded from the fi t due to the effect of CO in those 
bands. This supersedes the value of 18.39 ±0.06 found in lFreedman fc Mado rc (2010J) and is 



in excellent agreement with other independent measures rece ntly reviewed by IWalkerl (120111 ) 



who finds a composite distance modulus of 18.48 ± 0.05 and lLaney et al.l ( 120121 ) who find a 
red clump distance of 18.47 ± 0.02. 



7. Period-Color Relationships 

7.1. The Period-Color Diagram 

CO absorption in the 4.5 /iin band for cooler stars produces a significant period-color 
relation that moves the mean [3.6] — [4.5] color toward the blue at cooler temperatures, and 
longer periods. This effect is driven, as is the case for the light curve color variations (see 
below), by the temperature dependence of CO dissociation and not by the thermal color- 
temperature which has little effect on the slope of the continuum at these long wavelengths. 
The [3.6] — [4.5] Period-Color Relation is shown in Figure [6] and the weighted least squares 
fit to the de-reddened data (omitting Y Oph) is: 

[3.6] - [4.5] w = -0.09(±0.01)(logP - 1.0) - 0.03(±0.01). For comparison the LMC 



Period-Color Relation is [3.6]-[4.5] LM c = -0.09(±0.01)(logP-1.0)+0.01(±0.01) (IScowcroft et al 



201 lh . The slopes of these fits are consistent, but the zero-points differ by 0.04 ± 0.01 mag. 



Ngeow et al.l (120121 ) have presented a number of models for theoretical PC diagrams in 
the mid-IR. Their summary tables give slopes and zero-points for several models of the PLs 
and [3.6] — [4.5] color covering a range of Helium and metal abundance. Several of their 
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model PCs are plotted with the data in Figure [6) Comparison of the theoretical slope of the 
[3.6] — [4.5] PC diagram shows reasonably good agreement with the Galaxy color data for a 
Helium abundance of Y = 0.31. 



7.2. Color Curves and CO Absorption Models 



As pointed out in §4.11 the [3.6] — [4.5] color is characterized by systematic variations 
through the cycle. The color amplitude is also closely related to period, the longer period 
stars having the largest variations. This effect, in complete analogy with the cause of the 
Period-Color re l ation, is ag ain due to CO absorp tion in the 4.5 /im band, as discussed by 
Scowcroft et al.l ( 20111 ) and iMarengo et al.l ( 20101 ) . The color variation extends only toward 
bluer colors from a baseline red color limit of ~ 0.01 mag; this can be seen Figure [9j The blue 
extent (blue indicating more absorption at 4.5 /im) increases with period as the Cepheids 
reach intrinsically cooler temperatures. The effect of CO has only recently been observed 
over entire Cepheid pulsation cycles (see also IScowcroft et al.ll201ll ). 



To quantify the behavior of both the overall PC Relation and the color c urves, we have 



computed several synthetic spectra using appropriate Kurcz stellar models (IKuruczl Il993 



Castelli &: Kuruczll2003t ISbordone et al.l 12004 ; ISbordond 120051 ). Figure (TD] shows the results. 
They indicate that at temperatures greater than approximately 6000 K absorption due to 
CO is nearly non-existent. As the temperature falls below 6000 K CO absorption in the [4.5] 
band sets in, leading to the diminished flux observed in the 4.5 /im light curve. The result 
is that the color curves should have larger amplitudes for Cepheids with longer periods, as 
they reach intrinsically cooler temperatures. This is precisely the behavior exhibited in the 
observed color curves. For the shorter-period Cepheids the color amplitude is diminished 
because these Cepheids are intrinsically hotter and the CO remains dissociated over a longer 
portion of the pulsation cycle. 

Part of the systematic offset in the PC fits (Galaxy versus LMC) is plausibly explained 
by the difference in metal abundance between the two galaxies. Fig [TO] shows that this shift 
should amount to ~ 0.02 mag offset per 0.5 dex change in metal abundance. Other effects 
such as rotation may also play a role, and in any case we are dealing with a very small effect. 



8. Summary 



In this work we have presented the first results from the Galactic Cepheid campaign 
of the CHP. Light curves created from uniformly-spaced observations with high-precision 
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photometry yield intensity-mean magnitudes for 37 Galactic Cepheids spanning a range of 
periods from 4-70 days. 



Using the precise geometric parallax measurements from Benedict et al.l (120071 ) we have 
found a Galactic zero-point (set to log P=1.0) for the 3.6 /zm Period-Luminosity (Leavitt) 
Law of —5.80 ± 0.03 mag. Comparing this to the LMC zero-point we find an LMC distance 
modulus of 18.48 ± 0.04 mag, which is confirmed using a multi-wavelength analysis. The 
uncertainty represents a fac tor of 2 improvement over previous Key Project measurements 
(IFreedman fc Madorell2010l ) and will be made stronger with future geometric parallaxes to 
the full sample from GAIA. The impli cations of this revised L MC distance modulus on the 
Key Project distances are discussed in IFreedman et al.l (120121 ). 



The well-sampled light curves reveal a strong color variation for Cepheids with periods 
longer than 10 days. A second and related result is a clear Period-Color relation. Both 
correlations are caused by enhanced temperature-sensitive CO absorption at 4.67 /xm in 
longer period, intrinsically cooler Cepheids. 
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A. Corrections for Systematic Effects on the Photometry 

Saturated Pixels. The first step in the reduction is to find and mask markedly non-linear 
or saturated pixels. This is particularly important for stars as bright as the Cepheids in this 
program. Our routine works as follows: MOPEX fits the profile by masking unwanted pixels, 
and by assigning weights using the uncertainty image. Saturated pixels (near the center of 
the PSF) will tend to have lower dispersions and hence artificially high weights (in the limit 
of complete saturation the dispersion will be zero). We found valid upper thresholds by 
experiment: the dispersions for the nine dither positions of a saturated star were found as a 
function of threshold level and the best level chosen. The final upper threshold values were 
10,000 and 12,000 DN for the 3.6 and 4.5 /im bands, respectively. These were lower than 
than those recommended in the IRAC Handbook. 

Point Response Profiles and Pixel Phase Corrections. Each stellar profile was fit with 
a Point Response Function (PRF) profile, a procedure that minimizes the residual between 
the input frame and a standard PRF supplied by MOPEX. Rather than having a functional 
form, the PRF is a look-up table containing different representations of a point source at 
various pixel phases. These are the distances from the center of the stellar profile to the 
center of the nearest (integer valued) pixel. This complication arises because pixels do not 
have uniform response across them. The Pixel Phase Correction (PPC) was included in the 
PRF tables provided by the SSC for the cold mission. This was not the case for the warm 
mission. The correction for an arbitrary profile was thus determined using all the data to find 
a empirical PPC as follows. For every nine-point dither pattern constituting a measurement 
one has an average count, and nine deviations from that count. (The deviations arise from 
the pixel phase variation.) The left hand side of Figure [TT] shows all those deviations plotted 
as a function of pixel phase. The total number of points is 9 x 24 x 37 (9 dither positions per 
measurement, 24 light curve points per star, 37 stars). The strong correlation represents the 
residual PPC, which is easily removed to yield corrected data. The correction is of the form: 
fppc = u + v — p) , where p = a/ ((x — nint(x)) 2 + (y — nint(y)) 2 ) is the pixel phasdE 
The right-hand side of Figure [TTJ shows the results of applying the above correction to the 
data points on the left-hand side. The results are seen to be satisfactory. 

The PRF photometry methods correct for several other systematic effects: 

(1) The fits have slightly different coefficients for saturated and unsaturated data, a 
systematic effect that probably arises from centroid offsets in masking central pixels. The 



7 This equation is taken from the IRAC Handbook. The coefficients guarantee that the average deviation 
is zero and thus the correction does not introduce a spurious shift in the average measurement. 
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coefficients were confirmed by masking non-saturated data. 

(2) The PRF varies across the array, but a fixed PRF was used for the sub-array data, 
located at column 233, row 35. For the full frame data (CF Cas) a lookup table was used to 
find the nearest PRF at each position. 

(3) MOPEX reports fluxes at the center of pixel flux and is normalized to a radius of 10 
pixels, i.e., the IRAC standard aperture. Because the flux reported is that for a pixel phase 
of zero, an additional correction factor (f corr ) is applied to bring the flux to the average pixel 
phase, which corresponds to where the flux zero points were defined. 

Finally, the flux is placed on the standard Vega magnitude system by dividing by the 
photo metric flux zero-poi nts (zp); 280.9±4.1 and 179.7±2.6 [Jy], for [3.6] and [4.5], respec- 



tively ([Reach et al.ll2005l ). The final corrected magnitude (m) is found from the PRF flux 



(Fprf), which is reported in //Jy, by the following relation: 

m = —2.5 log 



Fprf ■ 10 6 



, fcorr ' fppc ' Zp 

Table [12] contains the constants for each channel. 

Image Persistence or Latency. The sub-array data were taken with short exposure 
times on bright objects, and with short settling times between dithers. As a consequence, 
the data are prone to short-term image persistence from previous dithers and observations; 
see Figure [12j A multi-stage process was undertaken to mitigate image persistence for each 
frame. First, the stellar profile was fit in each of the nine dithered frames using the PRF- 
fitting algorithm in the MOPEX script described above. The residual images were averaged 
together using a nearest neighbor weighting scheme for each dither position. Lowest weights 
were given to future frames and higher weights were given to the most recent frames to 
create an approximate map of the image persistence at each dither position. The persistence 
maps were then subtracted from the original data leaving only the source; see Figure [T3j 
Note that this process simply describes creating a local background frame by combining 
dithered frames where the source (s) has been modeled and subtracted rather than masked. 
Each background/persistence subtracted image was passed through the MOPEX pipeline to 
perform a second and final PRF fit. 

The latency can affect aperture photometry measurements by nearly five-percent in the 
short (0.02 s) sub-array data. The effect is larger for shorter exposure times when there is 
less time for the latency to dissipate. One advantage of the PRF-fitting algorithm is that it 
tends to ignore the latent pixels when fitting a PRF. After the latent image subtraction the 
PRF photometry changes by less than 1%. 

Photometry Checks. As a check of the photometric fidelity the standard star HD 165459 
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was processed in the manner described above for various exposure times in sub-array and 
full-array mode resulting in non-saturated and saturated data. To quantify the effect of 
persistence both aperture and PRF photometry were performed on the standard star both 
prior to, and after, the persistence correction. Prior to the correction, the aperture photom- 
etry was consistently reporting 5-15% higher flux than expected for the shortest exposure 
times while after the persistence correction the aperture photometry was typically only 0-5% 
higher than expected. The PRF photometry was relatively unaffected by the persistence cor- 
rection, changing by less than 1% before and after the persistence correction. The standard 
star comparisons represent the worse case scenario since the long exposures (saturated data) 
were taken just prior to the short exposures and subsequently suffered from a relatively large 
amount of persistence. The final PRF photometry results for the non-saturated standard 
data are 6.588 ± 0.007 mag and 6.571 ± 0.008 mag and 6.587 ± 0.014 mag and 6.564 ± 0.017 
mag for the saturated data. Both are in good agre ement with the stan dard magnitudes of 
6.593 ± 0.029 mag and 6.575 ± 0.028 mag found bv iReach et "aD J2005h . For the Cepheids, 
the difference in PRF photometry before and after the persistence correction was also less 
than 1%; except for i Car, the most severely saturated star in our sample, in which case the 
difference was 6% and 3% in Channels 1 and 2. We report here the final PRF photometry 
resulting from a persistence-subtracted image. 
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Table 1. CHP selected Galactic Cepheids and adopted parallaxes. 



ID 


log(P) 1 


R A 1 


DEC 1 


E(B-V) 2 


vr [no _l j 


7T 1 IVIO 1 


-TT-rTR^Rl 5 
vr [irxorjj 




[days] 




T9nnn 


[mags 






mac 


t > i .'I ic 

111(1.5 


S Vul 


1 


835 


19:48:23.8 


+27:17:11 


0.925 ± 


022 




0.32 ± 0.04 6 


0.27 + 0.01 


GY Sge 


1 


708 


19:35:13.6 


+19:12:09 


1.187 + 


170 






0.35 + 0.01 


SV Vul 


1 


653 


19:51:30.9 


+27:27:37 


0.508 + 


013 




0.56 + 0.03 


0.53 + 0.01 


U Car 


1 


589 


10:57:48.2 


-59:43:56 


0.265 + 


010 






0.71 ± 0.01 


I Car 


1 


551 


09:45:14.8 


-62:30:28 


0.154 + 


011 


2.01 + 0.20 




1.93 + 0.02 


T Mon 


1 


432 


06:25:13.0 


+07:05:09 


0.181 + 


010 






0.76 + 0.01 


WZ Sgr 


1 


339 


18:16:59.7 


-19:04:33 


0.501 + 


007 




0.56 + 0.03 


0.56 + 0.01 


RU Set 


1 


294 


18:41:56.4 


-04:06:38 


0.929 ± 


010 




0.60+0.03 


0.53 + 0.01 


SZ Aql 


1 


234 


19:04:39.4 


+01:18:22 


0.544 ± 


016 






0.47 + 0.01 


Y Oph 


1 


234 


17:52:38.7 


-06:08:37 


0.671 + 


008 






1.83 + 0.02 


CD Cyg 


1 


232 


20:04:26.6 


+34:06:44 


0.487 ± 


014 






0.41 ± 0.01 


X Cyg 


1 


214 


20:43:24.2 


+35:35:16 


0.241 + 


010 




0.82 ± 0.02 


0.89 + 0.01 


TT Aql 


1 


138 


19:08:13.7 


+01:17:55 


0.442 ± 


011 






1.03 + 0.02 


V340 Nor 


1 


053 


16:13:17.4 


-54:14:06 


0.313 + 


015 




0.58 + 0.04 


0.58 + 0.03 


TW Nor 


1 


033 


16:04:55.2 


-51:57:13 


1.172 + 


013 




0.52 ± 0.03 7 


0.46 + 0.02 


C Gem 


1 


007 


07:04:06.5 


+20:34:13 


0.018 + 


007 


2.78 + 0.18 


2.82 ± 0.12 8 


2.59 + 0.06 


Dor 





993 


05:33:37.5 


-62:29:23 


0.058 ± 


009 


3.14 + 0.16 




3.06 + 0.05 


S Nor 





989 


16:18:51.8 


-57:53:59 


0.177 + 


007 




1.11 ± 0.02 


1.05 + 0.01 


S Sge 





923 


19:56:01.3 


+16:38:05 


0.104 + 


008 






1.49 + 0.04 


DL Cas 





903 


00:29:58.6 


+60:12:43 


0.498 ± 


007 




0.60 + 0.01 




U Vul 





903 


19:36:37.7 


+20:19:59 


0.619 + 


009 






1.52 + 0.03 


W Sgr 





881 


18:05:01.2 


-29:34:48 


0.109 + 


007 


2.28 + 0.20 




4.61 ± 0.34 


r) Aql 





856 


19:52:28.4 


+01:00:20 


0.121 + 


008 






3.92 ± 0.06 


U Aql 





847 


19:29:21.4 


-07:02:39 


0.362 ± 


010 






1.69 + 0.05 


XSgr 





846 


17:47:33.6 


-27:49:51 


0.227 + 


013 


3.00 + 0.18 




3.11 ± 0.05 


U Sgr 





829 


18:31:53.3 


-19:07:30 


0.409 ± 


008 




1.73 + 0.08 


1.73 + 0.02 


V367 Set 





799 


18:33:35.2 


-10:25:38 


1.255 + 


016 




0.61 ± 0.01 




YSgr 





761 


18:21:23.0 


-18:51:36 


0.203 ± 


007 


2.13 + 0.29 




2.29 + 0.08 


V Cen 





740 


14:32:33.1 


-56:53:16 


0.282 ± 


008 




1.39 + 0.03 


1.65 + 0.12 


CV Mon 





731 


06:37:04.8 


+03:03:50 


0.733 ± 


014 




0.61 ± 0.02 


0.65 + 0.01 


8 Cep 





730 


22:29:10.3 


+58:24:55 


0.073 + 


007 


3.66 + 0.15 


3.61 + 0.20 9 


3.75 + 0.07 


CEa Cas 





711 


23:58:09.5 


+61:12:49 


0.549 ± 


010 




0.29 + 0.02 10 


0.32 + 0.02 11 


CF Cas 





688 


23:58:18.0 


+61:13:16 


0.549 ± 


010 




0.29 + 0.02 10 


0.32 + 0.02 11 


CEb Cas 





651 


23:58:09.1 


+61:12:49 


0.549 ± 


010 




0.29 + 0.02 10 


0.32 + 0.02 11 


FF Aql 





650 


18:58:14.7 


+ 17:21:39 


0.204 + 


008 


2.81 + 0.18 




2.71 ± 0.08 


T Vul 





647 


20:51:28.2 


+28:15:02 


0.073 + 


008 


1.90 + 0.23 




1.84 + 0.02 


RT Aur 





571 


06:28:34.1 


+30:29:35 


0.058 + 


011 


2.40 + 0.19 




2.57 + 0.03 



Note. — Three different methods of distance determination are examined: Hubble Space Telescope (HST) 
parallaxes, Main-Sequence (MS) fitting to candidate cluster and Infrared Surface Brightness (IRSB) method. 
The latter two methods were converted to parallax measures for comparison. 

1 Values adopted from the General Catalog of Variable Stars dSamus et al.ll2009h . The periods for V340 Nor 
and U Car were recomputed; see text. 

2 Values adopted from the average value in Table [2] 

3 HST parallaxes adopted from [Benedict et al. | [bOOTjl . The Lutz-Kelker bias correction is applied in the PL 
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analysis; see Table [5] 

4 Cluster parallaxes adopted from iTurnerl l l201Ci) . 
5 IRSB parallaxes adopted from lStorm et aL 
^Turner! j201lh . 
iMaiaess etal] d201lh . 
^Maiaess et al.l ll2012bh . 
^Maiaess et al ] ll2012ah . 
le lTurner fc Burke! fo00$) . 
^Matthews et alj l ll995l) . 
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Table 2. CHP selected Galactic Cepheids and reddenings. 



ID 


E(B-V)pft ot 


E(B- V) s pec 




E(B- V) s pace 




E(B-V)ai;e 


S Vul 





.727 ± 0.042 


0. 


940 ± 0.051 


1 


.020 ± 


0. 


030 





.925 ± 0.022 


GY Sge 


1 


.187 ± 0.170 














1 


.187 ± 0.170 


SV Vul 





.461 ± 0.022 


0. 


.510 ±0.020 





.590 ± 


0. 


.030 





.508 ± 0.013 


U Car 





.265 ± 0.010 

















.265 ± 0.010 


I Car 





.147 ±0.013 









.170 ± 


0. 


020 1 





.154 ± 0.011 


T Mon 





.181 ± 0.011 





.179 ±0.029 













.181 ± 0.010 


WZ Sgr 





.431 ± 0.011 


0. 


.458 ± 0.058 


o 


.560 ± 


0. 


.010 


o 


.501 ± 0.007 


RU Set 





.921 ± 0.012 









.950 ± 


0. 


.020 





.929 ± 0.010 


SZ Aql 





.537 ±0.017 


0. 


.588 ±0.041 













.544 ± 0.016 


Y Oph 





.645 ± 0.015 





.683 ±0.010 













.671 ± 0.008 


CD Cve 





.493 ±0.015 


0. 


.447 ± 0.040 










o 


.487 ± 0.014 


X Cvg 





.228 ±0.012 


0. 


.239 ±0.029 





.280 ± 


0. 


.020 


o 


.241 ± 0.010 


TT Aql 





.438 ±0.011 


0. 


.480 ± 0.036 













.442 ± 0.011 


V340 Nor 





.321 ± 0.018 









.290 ± 


0. 


.030 





.313 ± 0.015 


TW Nor 


1 


.157 ±0.014 






1 


.360 ± 


0. 


.050 


1 


.172 ± 0.013 


C Gem 





.014 ±0.011 





.031 ±0.041 





.020 ± 


0. 


010 1 





.018 ± 0.007 


Dor 





.052 ± 0.010 






o 


.080 ± 


0. 


020 1 


o 


.058 ± 0.009 


S Nor 





.179 ±0.009 


0. 


.268 





.170 ± 


0. 


.010 


o 


.177 ± 0.007 


S Sge 





.100 ±0.010 


0. 


.116 ±0.016 













.104 ± 0.008 


DL Cas 





.488 ± 0.010 


0. 


.487 ±0.024 


o 


.510 ± 


0. 


.010 


o 


.498 ± 0.007 


U Vul 





.603 ±0.011 


0. 


.663 ±0.018 













.619 ± 0.009 


W Sgr 





.108 ±0.011 


0. 


.079 ±0.017 





.120 ± 


0. 


010 1 





.109 ±0.007 


r\ Aql 





.130 ±0.009 





.096 ±0.015 













.121 ± 0.008 


U Aql 





.360 ±0.010 


0. 


.416 













.362 ±0.010 


XSgr 





.237 ±0.015 


0. 


.219 





.190 ± 


0. 


030 1 





.227 ±0.013 


U Sgr 





.403 ± 0.009 


0. 


.398 ± 0.022 





.500 ± 


0. 


.030 





.409 ± 0.008 


V367 Set 


1 


.231 ± 0.025 






1 


.270 ± 


0. 


.020 


1 


.255 ± 0.016 


YSgr 





.191 ± 0.010 


0. 


.182 ±0.021 





.220 ± 


0. 


010 1 





.203 ± 0.007 


V Cen 





.292 ±0.012 





.167 





.280 ± 


0. 


.010 





.282 ± 0.008 


CV Mon 





.722 ± 0.022 


0. 


.681 





.750 ± 


0. 


.020 





.733 ±0.014 


<5 Cep 





.075 ± 0.010 


0. 


.087 





.070 ± 


0. 


010 1 





.073 ± 0.007 


CEa Cas 






0. 


.503 













.549 ± 0.010 2 


CF Cas 





.553 ±0.011 


0. 


.527 ±0.025 













.549 ± 0.010 


CEb Cas 






0. 


.479 













.549 ± 0.010 2 


FF Aql 





.196 ±0.010 


0. 


.224 ±0.017 





.210 ± 


0. 


020 1 





.204 ± 0.008 


T Vul 





.064 ±0.011 


0. 


.068 ±0.015 





,110± 


0. 


020 1 





.073 ± 0.008 


RT Aur 





.059 ±0.013 


0. 


.050 ± 0.036 





.060 ± 


0. 


030 1 





.058 ±0.011 



Note. — Reddenings adopte d from the literature includ e: average photometric 
reddenings dFouque et al ]|2007l). s pectroscopic reddenin gs iKovtvukh et alll200Sf) 
and space reddenings llTurnerll201(it iBenedict et alj|2007h ; the last column contains 
the weighted average. When no error estimate is available for the spectroscopic 
reddening 0.05 is used. 

1 Space reddenings adopted from lBenedict et all ll2007l ) , otherwise iTurne'rl ll2010l) . 

2 CEa & CEb Cas are here assumed to have same excess as CF Cas. 
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Table 3. Spitzer IRAC photometry for Galactic Cepheids 



ID 


HMJD 1 


[3.6] 


[4.5] 




[days] 


[mag] 


[mag] 









1401 


■> 


■ Ulo 


EE 


n 

u 


UJ-U 


o 




EE 


n 

u, 




CD 


Cyg 


55147 


.8152 


5 


.659 


± 





.012 


5 


.740 


± 





.007 


CD 


Cyg 


55148 


.3628 


5 


.675 


± 





Oil 


5 


.728 


± 





.007 


CD 


Cyg 


55149. 


.3595 


5 


.723 


± 





.015 


5 


.730 


± 


0. 


.004 


CD 


Cyg 


55150. 


.2981 


5 


.713 


± 





.017 


5 


.696 


± 


0. 


.008 


CD 


Cyg 


55150. 


.8732 


5 


.587 


± 





.020 


5 


.555 


± 





.007 


CD 


Cyg 


55151 


.7063 


5 


.543 


± 





.018 


5 


.514 


± 


0. 


.008 


CD 


Cyg 


55152 


.0117 


5 


.528 


± 





.010 


5 


.505 


± 


0. 


.007 


CD 


Cyg 


55153 


.0145 


5 


.463 


± 





Oil 


5 


.446 


± 





.005 


CD 


Cyg 


55153 


.4658 


5 


.447 


± 





.016 


5 


.423 


± 


0. 


.008 


CD 


Cyg 


55154. 


.1197 


5 


.409 


± 





.007 


5 


.392 


± 


0. 


.006 


CD 


Cyg 


55154 


.7093 


5 


.381 


± 





.013 


5. 


.370 


± 





.008 


CD 


Cyg 


55155. 


.5184 


5 


.353 


± 





.015 


5 


.364 


± 





.007 


CD 


Cyg 


55156. 


.2294 


5 


.332 


± 





.007 


5 


.375 


± 


0. 


.004 


CD 


Cyg 


55157 


.2037 


5 


.328 


± 





.014 


5 


.392 


± 





007 


CD 


Cyg 


55157 


.8669 


5 


.331 


± 





.015 


5 


.412 


± 





008 


CD 


Cyg 


55158. 


.4195 


5 


.337 


± 





.015 


5 


.425 


± 


0. 


.005 


CD 


Cyg 


55158. 


.9929 


5 


.347 


± 





Oil 


5 


.450 


± 





.006 


CD 


Cyg 


55159 


.7584 


5 


.369 


± 





.007 


5 


.481 


± 





009 


CD 


Cyg 


55160 


.8400 


5 


.416 


± 





.010 


5 


.537 


± 


0. 


.008 


CD 


Cyg 


55161 


.1906 


5 


.430 


± 





.007 


5 


.562 


± 


0. 


.005 


CD 


Cyg 


55161 


.8165 


5 


.470 


± 





Oil 


5 


.602 


± 





.002 


CD 


Cyg 


55162. 


.8009 


5 


.529 


± 





.005 


5 


.662 


± 


0. 


.007 


CD 


Cyg 


55163 


.3453 


5 


.570 


± 





.007 


5 


.698 


± 


0. 


.004 



1 The heliocentric modified julian date (HMJD) was de- 
termined from averaging the 18 "HMJD" header keywords 
at each epoch (9 exposures in each Channel). Note: HMJD 
= HJD-2,400,000.5 

Note. — The quoted errors represent the random photo- 
metric errors determined from the variance of the dithered 
data. The systematic zero-point errors are 0.016 for both 
the [3.6] and [4.5] bands. The full table is available with the 
online version of this paper. 
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Table 4. Intensity averaged [3.6] and [4.5] apparent magnitudes for Galactic Cepheids 



ID 


loe-CPI 




^3 61 1 ) 1 




(\A 51 1 ) 1 


/[3 61-14 5}) 2 




\A p. v<^l 




L illtx 6°J 




[ryi a rrcl 


\vn Pi p^;1 
L iiltx 6 s J 


S Vul 


1.835 


4. 


.358 ± 0.003 


4. 


394 ± 0.003 


-0.036 ± 0.002 


GY Sge 


1.708 


4. 


.311 ± 0.003 


4. 


.345 ± 0.004 


-0.034 ± 0.003 


SV Vul 


1.653 


3 


.711 ± 0.005 


3 


.788 ±0.005 


-0.077 ±0.003 


U Car 


1.589 


3 


.357 ± 0.005 


3 


.415 ± 0.006 


-0.058 ± 0.003 


I Car 


1.551 





.925 ± 0.004 


1 


.047 ± 0.004 


-0.122 ± 0.003 


T Mon 


1.432 


3 


.359 ± 0.006 


3 


.425 ± 0.005 


-0.066 ± 0.003 


WZ Sgr 


1.339 


4. 


.364 ± 0.006 


4. 


.443 ± 0.005 


-0.079 ± 0.003 


RU Set 


1.294 


4. 


.856 ± 0.005 


4. 


.873 ± 0.005 


-0.017 ±0.002 


SZ Aql 


1.234 


4. 


.981 ± 0.005 


5 


.032 ± 0.005 


-0.051 ± 0.002 


Y Oph 


1.234 


2 


.528 ±0.002 


2 


.500 ± 0.002 


0.028 ±0.001 


CD Cyg 


1.232 


5. 


.477 ± 0.005 


5. 


.530 ± 0.005 


-0.053 ±0.002 


X Cyg 


1.214 


3 


.678 ± 0.005 


3 


.728 ±0.005 


-0.050 ± 0.002 


TT Aql 


1.138 


3 


.875 ± 0.005 


3 


.909 ± 0.004 


-0.034 ±0.002 


V340 Nor 


1.053 


5. 


.453 ± 0.002 


5 


.480 ± 0.002 


-0.027 ± 0.001 


TW Nor 


1.033 


6 


152 ± 0.003 


6 


160 ±0.004 


-0.008 ± 0.002 


C Gem 


1.007 


2 


025 ± 0.002 


2 


.037 ± 0.003 


-0.012 ±0.001 


fS Dor 


0.993 


1 


858 ± 0.003 


1 


.871 ± 0.003 


-0.013 ± 0.002 


S Nor 


0.989 


4. 


.066 ± 0.003 


4 


.085 ± 0.003 


-0.019 ± 0.002 


S Sge 


0.923 


3 


.652 ± 0.003 


3 


661 ± 0.003 


-0.009 ± 0.001 


DL Cas 


0.903 


5 


.783 ± 0.003 


5 


.791 ± 0.003 


-0.008 ± 0.001 


U Vul 


0.903 


3 


.797 ± 0.003 


3 


.778 ± 0.003 


0.019 ± 0.001 


W Sgr 


0.881 


2 


.721 ± 0.003 


2 


.719 ± 0.004 


0.002 ± 0.001 


r] Aql 


0.856 


1 


864 ± 0.003 


1 


.865 ± 0.003 


-0.001 ± 0.001 


U Aql 


0.847 


3 


.738 ± 0.003 


3 


.736 ± 0.003 


0.002 ± 0.001 


XSgr 


0.846 


2 


.423 ± 0.003 


2 


.409 ± 0.003 


0.014 ± 0.001 


U Sgr 


0.829 


3 


824 ± 0.003 


3 


.822 ± 0.003 


0.002 ± 0.001 


V367 Set 


0.799 


6 


.437 ± 0.001 


6 


.410 ± 0.002 


0.027 ± 0.001 


Y Sgr 


0.761 


3 


.486 ± 0.003 


3 


.483 ± 0.003 


0.003 ± 0.001 


V Cen 


0.740 


4. 


405 ± 0.003 


4. 


.400 ± 0.003 


0.005 ± 0.001 


CV Mon 


0.731 


6 


375 ± 0.003 


6 


.360 ± 0.003 


0.015 ± 0.001 


8 Ccp 


0.730 


2 


.221 ± 0.003 


2 


217 ± 0.003 


0.004 ± 0.001 


CEa Cas 


0.711 


7. 


.720 ± 0.003 


7. 


705 ± 0.003 


0.015 ± 0.001 


CF Cas 


0.688 


7. 


809 ± 0.002 


7. 


783 ± 0.002 


0.026 ± 0.001 


CEb Cas 


0.651 


7. 


.888 ± 0.003 


7. 


.851 ± 0.003 


0.037 ± 0.001 


FF Aql 


0.650 


3 


378 ± 0.001 


3 


353 ±0.001 


0.025 ± 0.001 


T Vul 


0.647 


4. 


114 ±0.003 


4. 


111 ± 0.003 


0.003 ± 0.001 


RT Aur 


0.571 


3 


.853 ± 0.003 


3 


.849 ± 0.003 


0.004 ± 0.001 



1 Shown with the expected random errors from the averaging algorithm; 
see text. The systematic errors are 0.016 mags for both the [3.6] and [4.5] 
data. 

2 Because the errors in Channel 1 and 2 arc correlated, the color error 
was calculated independently by using the error algorithm described in the 
text. 



Table 5. HST Parallax Cepheids With LKH 



Ccphcid 


logP 


m [3.6] 


m [4.5] 


LKH 


(m - M) 


^[3.6] 


A [4 .5] 


M [3 . 6] 




M [4 .5] 










I'm pi p n l 


[m pi p n l 


Pm pi p n l 


\rr\ pi p n l 


\vr\Fi p~l 


Ttiipi p~1 




Ivnfi p~1 




i Car 


1.551 


0.925 ± 0.004 


1.047 ± 0.004 


-0.08 


8.56 ±0.22 


0.031 ± 0.004 


0.024 ± 0.003 


-7.67 ± 0. 


22 


-7.54 ± 0. 


22 


C Gem 


1.007 


2.025 ± 0.002 


2.037 ± 0.003 


-0.03 


7.81 ± 0.14 


0.004 ± 0.001 


0.003 ± 0.001 


-5.79 ± 0. 


11 


-5.78 ± 0. 


11 


P Dor 


0.993 


1.858 ± 0.003 


1.871 ± 0.003 


-0.02 


7.54 ± 0.11 


0.012 ± 0.002 


0.009 ± 0.002 


-5.69 ± 0. 


11 


-5.67 ± 0. 


11 


W Sgr 


0.881 


2.721 ± 0.003 


2.719 ±0.004 


-0.06 


8.27 ± 0.19 


0.022 ± 0.003 


0.017 ± 0.002 


-5.57 ± 0. 


19 


-5.57 ± 0. 


19 


XSgr 


0.846 


2.423 ± 0.003 


2.409 ± 0.003 


-0.03 


7.64 ± 0.13 


0.046 ± 0.005 


0.035 ± 0.004 


-5.27 ± 0. 


13 


-5.27 ± 0. 


13 


YSgr 


0.761 


3.486 ± 0.003 


3.483 ±0.003 


-0.15 


8.51 ± 0.30 


0.041 ± 0.004 


0.032 ±0.003 


-5.06 ± 0. 


30 


-5.06 ± 0. 


30 


5 Cep 


0.730 


2.221 ± 0.003 


2.217 ± 0.003 


-0.01 


7.19 ±0.09 


0.015 ± 0.002 


0.011 ± 0.002 


-4.99 ± 0. 


09 


-4.99 ± 0. 


09 


FF Aql 


0.650 


3.378 ± 0.001 


3.353 ±0.001 


-0.03 


7.79 ±0.14 


0.041 ± 0.004 


0.032 ±0.003 


-4.45 ± 0. 


11 


-4.47 ± 0. 


11 


T Vul 


0.647 


4.114 ± 0.003 


4.111 ± 0.003 


-0.12 


8.73 ±0.26 


0.015 ± 0.002 


0.011 ± 0.002 


-4.63 ± 0. 


26 


-4.63 ± 0. 


26 


RT Aur 


0.571 


3.853 ± 0.003 


3.849 ± 0.003 


-0.05 


8.15 ±0.17 


0.012 ± 0.003 


0.009 ± 0.002 


-4.31 ± 0. 


17 


-4.31 ± 0. 


17 



Note. — The LKH factor has already been included the reported distance modulus, thus the absolute magnitude is: M = m — (m — M) — A. 
See T able [6] for distance moduli without LKH. The distance moduli for /3 Dor and W Sgr differ from the values tabulated in lBenedict et al] 
l|2007T ): the values reported here have been confirmed with Benedict (private communication). 



Table 6. HST Cepheids Without LKH 



Ccphcid 


logP 


"»[3.6] 


"l[4.5] 


LKH 


(to - M) 


^[3.6] 






M[ 3 .6] 






M [4 .5] 






[days] 


[mag] 


r i 
[magj 


[mag] 


[mag] 


[mag] 




r i 

[mag] 


[mag] 






[mag] 




I Car 


1.551 


0.925 ± 0.004 


1.047 ± 0.004 


0.00 


8.48 ± 0.22 


0.031 ± 0.004 





.024 ± 0.003 


-7.59 ±0. 


22 


-7 


.46 ±0. 


22 


£ Gem 


1.007 


2.025 ± 0.002 


2.037 ± 0.003 


0.00 


7.78 ± 0.14 


0.004 ± 0.001 





.003 ± 0.001 


-5.76 ±0. 


14 


-5 


.75 ±0. 


14 


fS Dor 


0.993 


1.858 ± 0.003 


1.871 ± 0.003 


0.00 


7.52 ± 0.11 


0.012 ± 0.002 





.009 ± 0.002 


-5.67 ± 0. 


.11 


-5 


.65 ± 0. 


11 


W Sgr 


0.881 


2.721 ± 0.003 


2.719 ±0.004 


0.00 


8.21 ± 0.19 


0.022 ± 0.003 





.017 ± 0.002 


-5.51 ± 0. 


19 


-5 


.51 ± 0. 


19 


XSgr 


0.846 


2.423 ± 0.003 


2.409 ± 0.003 


0.00 


7.61 ± 0.13 


0.046 ± 0.005 





.035 ± 0.004 


-5.24 ±0. 


13 


-5 


.24 ± 0. 


13 


YSgr 


0.761 


3.486 ± 0.003 


3.483 ±0.003 


0.00 


8.36 ± 0.30 


0.041 ± 0.004 





.032 ± 0.003 


-4.91 ±0. 


30 


-4 


.91 ±0. 


30 


8 Cep 


0.730 


2.221 ± 0.003 


2.217 ± 0.003 


0.00 


7.18 ± 0.09 


0.015 ± 0.002 





Oil ± 0.002 


-4.98 ±0. 


09 


-4 


.98 ±0. 


09 


FF Aql 


0.650 


3.378 ± 0.001 


3.353 ±0.001 


0.00 


7.76 ± 0.14 


0.041 ± 0.004 





.032 ± 0.003 


-4.42 ±0. 


14 


-4 


.44 ±0. 


14 


T Vul 


0.647 


4.114 ±0.003 


4.111 ± 0.003 


0.00 


8.61 ± 0.26 


0.015 ± 0.002 





Oil ± 0.002 


-4.51 ± 0. 


26 


-4 


.51 ± 0. 


26 


RT Aur 


0.571 


3.853 ± 0.003 


3.849 ± 0.003 


0.00 


8.10± 0.17 


0.012 ± 0.003 





.009 ± 0.002 


-4.26 ±0. 


17 


-4 


.26 ±0. 


17 



Table 7. Cluster Cepheids 



Cepheid 


logP 


m [3.6] 


™[4.5] 


(m - M) 




A [3.6] 


^[4.5] 


A^[3.6] 




M [4 .5] 




[days] 


[mag] 


[mag] 


[mag] 






[mag] 


[mag] 


[mag] 




[mag] 


O VUl 


1 QOC 
l.OOO 


4.O0O ± U.UUo 


a QO/i J_ n nnQ 
4.oy4 it U.UUo 


12.47 ± 


.27 


n 
U. 


1 qq _!_ n m n 

loo it u.uiy 


n i/i/i J_ n m t; 
U. 144 it U.Ulo 


q Qn _i_ n 
— o.oU it U. 


OT 7 
Z I 


Q OO J_ n OT/ 

— o.ZZ ± V.Z t 


O V VUl 


1 fit: O 

l.OOo 


o 711 _i_ n nnt; 
O. / 11 it u.uuo 


q 7QO _i_ n nnt; 
o. loo it U.UUo 


11.26 ±0. 


.12 


n 
U. 


1 HQ J. n nil 
lUo it U.U11 


n H7Q 4- n nnc 
u.u* y it u.uuo 


i «t; _i_ n 
— i .00 it u. 


1 
. LZ 


1 t; t; _i_ n io 
— f .00 it U. 1Z 


WZj ogr 


1 QQO 

l.ooy 


4.oo4 ± U.UUo 


a 4 41 I n nnc; 
4.44o it U.UUO 


11.26 ± 0. 


.12 


U. 


1 no J_ n m n 
WZ it U.U1U 


n (\"7Q _i_ n nno 
U.U/o it U.UUo 


i nn _l n 
— f .UU it u. 


1 O 
. LZ 


a co j_ n i o 
— o.oy it U. vz 


T3TT Q,"f 


i.zy4 


/i Qtrfi _i_ n nnt; 
4.00D it U.UUO 


a av! _i_ n nne; 

4.o I o it U.UUO 


11.11 ± 0. 


.11 


n 
u. 


i qo J. n mo 
loy it u.uiy 


n 1K 4-n mi; 
U. 140 it U.Ulo 


— 0.44 it U. 


1 1 
. 1 1 


r qq _i_ n 1 1 
— o.oo it U. 11 




1 Q1/I 

1.Z14 


o.O / o ± U.UUo 


q 7oo j_ n nnt; 
o. / Zo ± U.UUo 


10.43 ± 0. 


.05 


U. 


c\a o _t_ n nnt; 
U4y it U.UUO 


n nQC _i_ n nn/i 
U.Uoo it U.UU4 


on J_ n 
— O.oU it u. 


nt; 
Uo 


a i a _i_ n nt; 
— 0. / 4 ± U.Uo 


V o4U IN or 


l.UOo 


c a co _i_ n nno 
0.40O it U.UU.Z 


t; a on _i_ n nno 

0.4oU it U.UU.Z 


11.18 ±0. 


.15 


U. 


nfi/i -t n nnv 
UD4 it U.UU / 


n n/io -t n nnt; 
U.U4y it U.UUO 


t; 70 -t n 


i t; 

10 


t; 7t; j. n it; 
— 0. / it U. 10 


TW Nor 


1.033 


6.152 ± 0.003 


6.160 ± 0.004 


11.42 ±0. 


.13 


0. 


238 ± 0.024 


0.183 ± 0.018 


-5.51 ± 0. 


.13 


-5.44 ±0.13 


£ Gem 


1.007 


2.025 ± 0.002 


2.037 ± 0.003 


7.75 ±0. 


.09 


0. 


004 ± 0.001 


0.003 ± 0.001 


-5.73 ±0. 


09 


-5.71 ± 0.09 


S Nor 


0.989 


4.066 ± 0.003 


4.085 ± 0.003 


9.77 ± 


,04 


0. 


036 ± 0.004 


0.028 ± 0.003 


-5.74 ± 0. 


,04 


-5.72 ±0.04 


DL Cas 


0.903 


5.783 ± 0.003 


5.791 ± 0.003 


11.11 ± 


.04 





101 ± 0.010 


0.078 ± 0.008 


-5.43 ±0. 


04 


-5.40 ± 0.04 


U Sgr 


0.829 


3.824 ± 0.003 


3.822 ± 0.003 


8.81 ± 


10 


0. 


083 ± 0.008 


0.064 ± 0.007 


-5.07 ± 0. 


.10 


-5.05 ±0.10 


V367 Set 


0.799 


6.437 ±0.001 


6.410 ± 0.002 


11.07 ± 0. 


.04 


0. 


255 ± 0.026 


0.196 ± 0.020 


-4.89 ±0. 


.04 


-4.86 ±0.04 


V Cen 


0.740 


4.405 ± 0.003 


4.400 ± 0.003 


9.28 ±0. 


.05 


0. 


057 ± 0.006 


0.044 ± 0.005 


-4.94 ± 0. 


.05 


-4.93 ±0.05 


CV Mon 


0.731 


6.375 ± 0.003 


6.360 ± 0.003 


11.07 ± 0. 


.07 


0. 


149 ±0.015 


0.114 ±0.012 


-4.85 ±0. 


.07 


-4.83 ±0.07 


<5 Cep 


0.730 


2.221 ± 0.003 


2.217 ± 0.003 


7.21 ± 0. 


.12 


0. 


015 ± 0.002 


0.011 ± 0.002 


-5.01 ± 0. 


.12 


-5.01 ± 0.12 


CEa Cas 


0.711 


7.720 ± 0.003 


7.705 ± 0.003 


12.69 ±0 


15 


0. 


111 ± 0.011 


0.086 ± 0.009 


-5.08 ±0. 


,15 


-5.07 ± 0.15 


CF Cas 


0.688 


7.809 ± 0.002 


7.783 ± 0.002 


12.69 ±0. 


.15 


0. 


111 ± 0.011 


0.086 ± 0.009 


-4.99 ±0. 


.15 


-4.99 ±0.15 


CEb Cas 


0.651 


7.888 ± 0.003 


7.851 ± 0.003 


12.69 ±0 


.15 


0. 


111 ± 0.011 


0.086 ± 0.009 


-4.91 ±0. 


.15 


-4.92 ±0.15 



Table 8. IRSB Cepheids 



Cepheid 


logP 


m [3.6] 


"l[4.5] 


(m - M) 




A [3.6] 




^[4.5] 


M [3 .6] 


M [4 .5] 




[days] 


[mag] 


[mag] 


[mag] 




[mag] 




[mag] 


[mag] 


[mag] 


O VUl 


l.ooo 


4.o0o ± U.UUo 


a QO/i J_ n nnQ 
4.oy4 it U.UUo 


IO Q A _!_ n HQ 

1Z.04 it U.Uo 


n 
U. 


1 oq j_ n m 
. loo it u.uiy 


n 
U. 


1 4 \ in m tr 
.144 it U.U10 


q pi J_ n no 
— 0.0/ it U.Uo 


q c;n _l n no 

—o.oy ± u.uo 


Vjj I ogc 


1 70.Q 
1. I UO 


A Q1 1 4-0. nnQ 
^i.Oll it u.uuo 


A oa k _i_ n r\r\A 
^.O^tO it u.uu^t 


1 os 4- n or 

IZ.Zo It U.UO 


u. 


OA1 4- flzlO 
Z^l It u.u^tz 


n 
U. 


1 4- flQO 
100 it u.uoz 


Q 01 4- fl HQ 
— O.Z1 It U.UO 


Qi9j.f1 07 
— 0. 1Z I U.Ul 


O V VUl 


1 ftt; q 

l.ooo 


q 7i 1 j_ n nnt; 
o. ( 11 ± U.UUO 


q 700 _i_ n nnt; 
0. loo ± U.UUO 


1 1 qq j_ n n/i 
11.00 it U.U4 


n 

U. 


1 nQ J_ n m 1 
lUo it U.U11 


n 
U. 


nvo J_ n nnQ 

u / y ± u.uuo 


7 77 -J- n f\A 

— ( . 1 ( it U.U4 


7 ft7 J_ n (\A 

— ( .0 / it U.U4 


TT Per 

u oar 


1 t;QO 

l.ooy 


q QK7 -I- n nnt; 

O.oO i it U.UUO 


q a 1 c _i_ n nnR 
O.410 it U.UUO 


1 n 7/i 4- n nQ 
1U. 1 4 it U.Uo 


u. 


n^A 4- n nnfi 
U04 it U.UUO 


U. 


C\A 1 4- n C\f\A 
U41 it U.UU4 


7 AA 4- n nQ 
— / .44 it U.Uo 


7 Q7 4- n nQ 
— / .0 / it U.Uo 


t oar 


1 EC 1 
1.001 


o not; j_ n ori/1 
U.yzo ± U.UU4 


1 n/i 7 _l n nn/i 
1.U4 f it U.UU4 


q c7 j_ n no 
o.Of it U.UZ 


0. 


031 ± 0.004 


n 

U. 


noA -L n nnQ 
Uz4 it U.UUo 


7 pq -L n no 
— / .00 it U.Uz 


7 t; c; j_ n no 
— / .00 it v.UZ 


T Mon 


1 /I QO 


o.ooy ± U.UUO 


q a oc; j_ n nnt; 
O.4Z0 it U.UUO 


1 n fin J_ n nQ 
1U.0U it U.Uo 





037 ± 0.004 


n 

U. 


noQ J_ n nnQ 
Uzo it U.UUo 


7 07 J_ n nQ 
— i .Zi it U.Uo 


1 on J_ n nQ 
— / .ZD ± U.Uo 


XKT7 Cerr- 

vv Zi ogr 


l.ooy 


zi Qfi/i _i_ n nnfi 
4.0D4 it U.UUO 


a a a q 4- n nnt; 

4.44o it U.UUO 


1 1 ofi 4- n n/i 
11. ZD it U.U4 


0. 


102 ± 0.010 


U. 


n7Q 4- n nnQ 
U10 ± U.UUo 


7 nn 4- n r\A 
— 1 .UU it U.U4 


ft qo 4- n n/i 
— O.oy it U.U4 


rt U OCT 


i.zy4 


a qkp _i_ n nnt; 

4.O00 it U.UUO 


a q7q _i_ n nnt; 
4.0 1 it U.UUO 


1 1 qo _i_ n n/i 
11.00 it U.U4 


0. 


189 ± 0.019 


U. 


1 /it; _l n m t; 
140 it U.U10 


fi 71 4- n nt; 
— O./l it U.UO 


ft fit; _i_ n n/i 
— 0.00 it U.U4 


C7 A 

oZi J\C\Y 


1 OQ/1 
1.Z04 


zi oqi 4- n nnt; 
4.yei it u.uuo 


c; nQO 4- n nnt; 
O.Uoz it U.UUO 


1 1 Q.A _i_ n nt; 
11.04 it U.UO 





110 ±0.012 


n 
u. 


not; _i_ n nno 
uoo it u.uuy 


ft 77 4- n nt; 
— O. ( ( ± U.UO 


ft fto 4- n nt; 
— o.oy it u.uo 




1 OQ/1 
1.Z04 


c:oo _i_ n nno 
Z.Ozo it U.UUz 


cenn 4- n nno 
Z.OUU it u.uuz 


q AO 4- n no 
o.oy it u.uz 


0. 


136 ±0.014 


n 

U. 


1 nc; 4- n m 1 
1U0 it U.Ull 


ft Qn 4- n nQ 
— O.oU it U.Uo 


ft 00 4- n nQ 
— O.zy it U.Uo 


OJJ oyg 


l.zoz 


t; AT7 J- n nnt; 
0.4 ( ( ± U.UUO 


c; kqh J_ n nnt; 
O.OoU it U.UUO 


1 1 cia J_ n nc; 
11. y4 it U.UO 


0. 


099 ±0.010 


n 

U. 


c\"7P J_ n nnQ 
U ( it U.UUo 


ft c;c _i_ n nt; 
— 0.00 it U.UO 


ft a q J_ n nt; 
— 0.4o it U.UO 


A ^yg 


1 01 /I 


q fi7o _i_ n nnt; 
0.0 1 It U.UUO 


700 _i_ n nnt; 

0. / Zo It U.UUO 


1 n oc; 4- n no 
1U.ZO it u.uz 





049 ± 0.005 


n 

u. 


nQQ 4- n nn/i 

Uoo It U.UU4 


ft fto 4- n nQ 
— O.Oz it U.Uv5 


ft rzp _i_ n nQ 
— 0.00 it U.Uo 


1 1 Aql 


1. loo 


q 07c: _t_ n one; 
0.0 / it U.UUO 


q ono J_ n nn/i 

o.yuy it u.uu4 


CkA j_ n n/i 

y.y4 ± U.U4 





090 ± 0.009 


n 
U. 


nfio J_ n nn7 
uoy it U.UU / 


ft 1 t; _L n C\A 
— 0.10 it U.U4 


ft 1 n j_ n (\a 
— 0. 1U it U.U4 


V o4U IN or 


1 nt;Q 

1.U00 


0.40o it U.UUz 


k Ann 4- n nno 

0.4oU it U.UUZ 


11 iQ4_nii 

11. lo It U. 1 1 


0. 


064 ± 0.007 


U. 


n/io 4- n nnt; 
U4y it u.uuo 


c; 70 4_ n 1 1 
— 0. / y it U.ll 


t; 7t; _i_ n 1 1 
— 0. (0 ± U.ll 


i vv in or 


l.UOO 


ft 1 ^9 4- n nnq 

O.IOZ It U.UUO 


« 1 «n _i_ n nn/i 
O.lOU it u.uu^i 


1 1 fiQ _i_ n no 
11. oy it u.uy 


0. 


.238 ± 0.024 


n 
U. 


1 qq 4- n m Q. 
100 It U.UIO 


c; 77 ±n in 

— O.I 1 It U. 1U 


t; 71 _|_ n 1 n 
— O. 1 1 it U. 1U 




i nn7 

l.UU f 


not; _l n nno 
z.Uzo it U.UUz 


nQ7 4- n nnQ 
.Z.Uo / it U.UUo 


7 oq 4- n nt; 
1 . y it u.uo 


0.004 ±0.001 


n 

U. 


nnQ 4- n nm 
UUo it U.UUl 


t; 01 4- n nt; 

— o.yi it u.uo 


t; on 4- n nt; 
— o.yu it U.UO 


p uor 


u.yyo 


1 Qcro 1 n nnQ 
I.0O0 it U.UUo 


1 071 _i_ n nnQ 
l.oi 1 it U.UUo 


7 t^7 4- n CiA 
/ .0 1 it U.U4 


0. 


012 ± 0.002 


U. 


nno -i- n nno 
uuy it U.UUz 


t; 7q _i_ n c\a 
— 0. / it U.U4 


c 71 4. n n/i 

— 0.(1 it U.U4 


o in or 


u.yoy 


A nftft 4- n nnQ 
4.U00 it U.UUo 


a net; _i_ n nnQ 
4.UoO it U.UUo 


qo 4- n no 
y.oy it v.yjz 


0. 


036 ± 0.004 


n 

U. 


noc 4- n nnQ 
Uzo it U.UUo 


t; Qfi 4- n no 
— O.oO it U.Uz 


c; Q/i 4- n no 

— O.o4 it V.uZ 


o bge 


u.yzo 


q ftt;o J_ n nnQ 
o.OOZ it U.UUo 


q fifii j_ n nnQ 
o.OOl it U.UUo 


1 q j_ n nfi 
y. lo it U.UO 


0. 


021 ± 0.003 


n 
U. 


m p J_ n nno 
U10 it U.UUz 


t; c;n J_ n nft 
— O.OU it U.UO 


<z ac\ i n nft 

— o.4y ± u.uo 


TT 

U VUl 


n ooq 

u.yuo 


q 707 4- n nnQ 
0. /y < it U.UUo 


q 770 _i_ n nnQ 
0.1 ( it U.UUo 


no 4- n c\A 

y.uy it U.U4 


0. 


126 ±0.013 


U. 


no7 4- n ni n 
uy ( it u.uiu 


£ ao 4- n n/i 

— 0.4Z it U.U4 


k /M 4_ n n/i 
— 0.41 it U.U4 


W ogr 


U.ool 


701 J_ n nnQ 
z,/zl it U.UUo 


71 J_ n nn/i 

z. ( iy ± U.UU4 


O.Oo it U. 10 


0. 


022 ± 0.003 


n 
U. 


m 7 j_ n nno 
Ul / it U.UUz 


q oq 1 n 1 ft 
— o.yo it U. 10 


Q OQ 1 n 1 ft 

— o.yo it U. 10 


fj Aql 


n Qt;ft 
U.oOO 


1 qpa 4- n nnQ 
1.504 it U.UUo 


1 QiRt; _i_ n nnQ 
l.oOO it U.UUo 


7 nQ 4- n nQ 
/ .Uo it U.Uo 


0. 


.025 ± 0.003 


n 

U. 


m 4- n nno 
uiy it U.UUZ 


k 1 4- n nQ 
— o.iy it u.uo 


; in ± n nQ 

— 0. iy it U.Uo 


U Aql 


0.847 


too t n nno 

3.738 ± U.UUo 


70/ 1 1 n nno 

3.736 ± 0.003 


or* 1 n dc 

8.86 ± 0.06 


0. 


073 ± 0.008 


0. 


056 ± 0.006 


tr nn I n n/"* 

— 5.20 ± 0.06 


— 5.18 ± 0.06 


XSgr 


0.846 


2.423 ±0.003 


2.409 ± 0.003 


7.54 ± 0.03 


0. 


046 ± 0.005 


0. 


035 ± 0.004 


-5.16 ± 0.04 


-5.16 ±0.04 


U Sgr 


0.829 


3.824 ±0.003 


3.822 ± 0.003 


8.81 ± 0.03 





.083 ± 0.008 


0, 


064 ± 0.007 


-5.07 ± 0.03 


-5.05 ±0.03 


YSgr 


0.761 


3.486 ±0.003 


3.483 ± 0.003 


8.20 ± 0.08 





041 ± 0.004 


0, 


032 ± 0.003 


-4.76 ± 0.08 


-4.75 ±0.08 


V Cen 


0.740 


4.405 ± 0.0U3 


4.400 ± 0.003 


8.91 ± 0.16 


0. 


057 ± 0.006 


0. 


044 ± 0.005 


-4.56 ± 0.16 


-4.56 ±0.16 


CV Mon 


0.731 


6.375 ±0.003 


6.360 ± 0.003 


10.94 ±0.03 


0. 


149 ±0.015 


0. 


114 ±0.012 


-4.71 ± 0.04 


-4.69 ±0.04 


5 Cep 


0.730 


2.221 ± U.0U3 


2.217 ± 0.003 


7.13 ±0.04 


0. 


.015 ±0.002 


0. 


011 ± 0.002 


-4.92 ± 0.04 


-4.92 ±0.04 


CEa Cas 


0.711 


7.72U±U.0U3 


7.705 ± 0.003 


12.47 ± 0.14 


0. 


111 ± 0.011 


0. 


086 ± 0.009 


-4.87 ± 0.14 


-4.85 ±0.14 


CF Cas 


0.688 


7.809 ± 0.002 


7.783 ± 0.002 


12.47 ± 0.14 


0. 


111 ± 0.011 


0, 


086 ± 0.009 


-4.78 ± 0.14 


-4.78 ±0.14 


CEb Cas 


0.651 


7.888 ±0.003 


7.851 ± 0.003 


12.47 ± 0.14 





111 ± 0.011 


0, 


086 ± 0.009 


-4.70 ± 0.14 


-4.71 ± 0.14 


FF Aql 


0.650 


3.378 ±0.001 


3.353 ± 0.001 


7.84 ± 0.06 


0. 


041 ± 0.004 


0. 


032 ± 0.003 


-4.50 ± 0.06 


-4.51 ± 0.06 


T Vul 


0.647 


4.114 ±0.003 


4.111 ± 0.003 


8.68 ± 0.02 


0. 


015 ± 0.002 


0. 


011 ± 0.002 


-4.58 ±0.02 


-4.58 ±0.02 



Table 8 — Continued 



Ccphcid log P 


m [3.6] 


"l[4.5] 


(m - M) 


A [3.6\ 


A [4 .5] 


M[3.6] 


^[4.5] 


[days] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 



RT Aur 0.571 



3.853 ±0.003 3.849 ± 0.003 7.95 ± 0.03 0.012 ± 0.003 0.009 ± 0.002 



-4.11 ± 0.03 



-4.11 ± 0.03 



Table 9. Period-Luminosity (Leavitt) Laws 



















wei 


rhts= -i 
0.1 2 






wei 


ghts = 

CT njr 

M 






weights 


1 






Galaxy 




Sample 






Band 


N 


(slope) a 


(intercept) b 




(T 


(slope) a 


(intercept) b 




u 


(slope) a 


(intercept) b 




a 


MW 


"HST 








[3 


6] 


10 


-3.40 ± 0.12 


-5.81 ± 0.04 





09 


-3.35 ± 0.22 


-5.80 ± 0.06 





10 


-3.39 ± 0.26 


-5 


80 ± 0.07 





09 


MW 


"HST 


6 < P < 60 






[3 


6] 


5 


-3.33 ± 0.18 


-5.82 ± 0.05 





09 


-3.33 ± 0.35 


-5.79 ± 0.06 





08 


-3.34 ± 0.39 


-5 


80 ± 0.08 





08 


MW 


"HST 


no LKH 






[3 


6] 


10 


-3.40 ± 0.12 


-5.75 ± 0.04 





09 


-3.33 ± 0.22 


-5.76 ± 0.06 





10 


-3.37 ± 0.26 


-5 


76 ± 0.07 





10 


MW 


^HST 


no LKH 6 < P 


< 


(ill 


[3 


6] 


5 


-3.27 ± 0.18 


-5.78 ± 0.05 





07 


-3.27 ± 0.35 


-5.75 ± 0.06 





06 


-3.28 ± 0.39 


-5 


76 ± 0.08 





07 


MW 


"MS 








[3 


6] 


18 


-3.00 ± 0.07 


-5.76 ± 0.02 





19 


-3.43 ± 0.08 


-5.77 ± 0.02 





17 


-3.11 ± 0.12 


-5 


75 ± 0.03 





19 


MW 


"MS 6 


< P < 60 






[3 


6] 


11 


-3.30 ± 0.12 


-5.68 ± 0.03 





19 


-3.70 ± 0.10 


-5.74 ± 0.02 





17 


-3.37 ± 0.17 


-5 


70 ±0.04 





19 


MW 


"IRSB 








[3 


6] 


34 


-3.42 ± 0.05 


-5.73 ± 0.02 





1.3 


-3.39 ± 0.02 


-5.74 ± 0.01 





13 


-3.41 ± 0.06 


-5 


73 ± 0.02 





13 


MW 


"IRSB 


6 < P < 60 






[3 


6] 


23 


-3.35 ± 0.08 


-5.75 ± 0.02 





1.3 


-3.31 ± 0.03 


-5.77 ± 0.01 





14 


-3.33 ± 0.09 


-5 


75 ± 0.03 





13 


LMC 1 


6 < P 


< 60 






[3 


6] 


80 


-3.31 ± 0.05 


12.70 ± 0.02 





11 


-3.37 ± 0.07 


12.71 ± 0.02 





11 


-3.31 ± 0.05 


12. 


70 ± 0.02 





11 


MW 2 


"HST 








[3 


6] 


10 




-5.80 ± 0.03 





10 




-5.79 ± 0.05 





10 




-5 


79 ± 0.06 





10 


MW 2 


"HST 


6 < P < 60 






[3 


6] 


5 




-5.82 ± 0.04 





09 




-5.79 ± 0.06 





08 




-5 


80 ± 0.08 





08 


MW 2 


^HST 


no LKH 






[3 


6] 


10 




-5.74 ± 0.03 





10 




-5.76 ± 0.05 





10 




-5 


75 ± 0.06 





10 


MW 2 


■""HST 


no LKH 6 < P 


< 60 


[3 


6] 


5 




-5.77 ± 0.04 





07 




-5.75 ± 0.06 





06 




-5 


76 ± 0.08 





07 


MW 2 


X MS 








[3 


6] 


18 




-5.75 ± 0.02 





21 




-5.76 ± 0.02 





17 




-5 


76 ± 0.03 





19 


MW 2 


"MS 6 


< P < 60 






[3 


6] 


11 




-5.67 ± 0.03 





19 




-5.74 ± 0.02 





18 




-5 


71 ± 0.04 





19 


MW 2 


""IRSB 








[3 


6] 


34 




-5.74 ± 0.02 





14 




-5.75 ± 0.01 





13 




-5 


74 ± 0.02 





13 


MW 2 


"IRSB 


6 < P < 60 






[3 


6] 


23 




-5.75 ± 0.02 





13 




-5.77 ± 0.01 





14 




-5 


75 ± 0.02 





13 


MW 


"HST 








[1 


5] 


10 


-3.26 ± 0.12 


-5.78 ± 0.04 


(1 


09 


-3.23 ± 0.22 


-5.77 ± 0.06 





09 


-3.26 ± 0.26 


-5 


77 ± 0.07 





09 


MW 


"HST 


6 < P < 60 






[4 


■ r >] 


5 


-3.13 ± 0.18 


-5.80 ± 0.05 





08 


-3.15 ± 0.35 


-5.77 ± 0.06 





08 


-3.15 ± 0.39 


-5 


77 ± 0.08 





08 


MW 


""HST 


no LKH 






H 


■ r >l 


10 


-3.27 ± 0.12 


-5.72 ± 0.04 





09 


-3.21 ± 0.22 


-5.74 ± 0.06 





10 


-3.25 ± 0.26 


-5 


73 ± 0.07 





10 


MW 


""HST 


no LKH 6 < P 


< 


60 


[4 


■ r >) 


5 


-3.07 ± 0.18 


-5.76 ± 0.05 





07 


-3.09 ± 0.35 


-5.73 ± 0.06 





06 


-3.09 ± 0.39 


-5 


74 ± 0.08 





06 


MW 


"MS 








[4 


■ r >) 


18 


-2.91 ± 0.07 


-5.72 ± 0.02 





19 


-3.33 ± 0.08 


-5.73 ± 0.02 





17 


-3.01 ± 0.12 


-5 


71 ± 0.03 





19 


MW 


IMS 6 


< P < 60 






[4 


■ r >] 


11 


-3.19 ± 0.12 


-5.64 ± 0.03 





19 


-3.61 ± 0.10 


-5.71 ± 0.02 





17 


-3.27 ± 0.17 


-5 


66 ± 0.04 





20 


MW 


"IRSB 








[4 


5] 


34 


-3.32 ± 0.05 


-5.70 ± 0.02 





1.3 


-3.28 ± 0.02 


-5.71 ± 0.01 





13 


-3.31 ± 0.06 


-5 


70 ± 0.02 





13 


MW 


"IRSB 


6 < P < 60 






[4 


■ r >) 


23 


-3.23 ± 0.08 


-5.72 ± 0.02 





12 


-3.18 ± 0.03 


-5.74 ± 0.01 





12 


-3.22 ± 0.09 


-5 


72 ± 0.03 





12 


LMC 1 


6 < P 


< 60 






[4 


■ r >l 


80 


-3.21 ± 0.06 


12.69 ± 0.02 





12 


-3.28 ± 0.08 


12.71 ± 0.03 





12 


-3.21 ± 0.06 


12 


70 ± 0.02 





12 


MW 2 


"HST 








[4 


5] 


10 




-5.77 ± 0.03 





10 




-5.77 ± 0.05 





09 




-5 


76 ± 0.06 





10 


MW 2 


""HST 


6 < P < 60 






[4 


■ r >l 


5 




-5.79 ± 0.04 





09 




-5.77 ± 0.06 





08 




-5 


77 ± 0.08 





08 


MW 2 


"HST 


no LKH 






[4 


5] 


10 




-5.71 ± 0.03 





10 




-5.73 ± 0.05 





10 




-5 


73 ± 0.06 





10 


MW 2 


"HST 


no LKH 6 < P 


< 60 


[4 


5] 


5 




-5.75 ± 0.04 





08 




-5.73 ± 0.06 





07 




-5 


74 ± 0.08 





07 


MW 2 


"MS 








[4 


5] 


18 




-5.72 ± 0.02 





21 




-5.72 ± 0.02 





18 




-5 


72 ± 0.03 





20 


MW 2 


7T MS 6 


< P < 60 






[4 


5] 


11 




-5.63 ± 0.03 





19 




-5.70 ± 0.02 





18 




-5 


67 ± 0.04 





20 


MW 2 


•""IRSB 








[4 


■ r >] 


34 




-5.71 ± 0.02 





1.3 




-5.72 ± 0.01 





13 




-5 


71 ± 0.02 





13 


MW 2 


"IRSB 


6 < P < 60 






[4 


■ r >] 


23 




-5.72 ± 0.02 





12 




-5.73 ± 0.01 





13 






72 ± 0.02 





12 



"'"The LMC sample does not contain Cepheids with periods less than 6 days. The LMC data are discussed in Scowcroft ct al. (2011). 
2 Force fit the LMC slope from the un-wcighting method to find the zero-point. 



Note. — The form of the PL relation used in these fits is: M — a{log P — 1) + b. The values in bold indicate our adopted PL slope and zero-points. To eliminate asymmetric 
rounding errors when reporting two significant figures in the zero-points a floor rounding function was used which rounds toward negative infinity rather than away from zero. 



Table 10. Multi- Wavelength parameters for Galactic Cepheids. 



Cepheid 


logP 


E(B-V) 


7r [mas] 


(m-M) 


LKH 


rriB 


mv 


171/ 


mj 


m H 


m K 


m [3.6] 


™[4.5] 


t Car 


1.551 


0.154 ±0.011 


2.01 ± 0.20 


8.48 ± 0.22 


-0.08 


4.986 


3.722 


2.552 


1.674 


1.198 


1.076 


0.925 


1.047 


C Gem 


1.007 


0.018 ±0.007 


2.78 ±0.18 


7.78 ± 0.14 


-0.03 


4.701 


3.889 


3.107 


2.585 




2.114 


2.025 


2.037 


P Dor 


0.993 


0.058 ± 0.009 


3.14 ±0.16 


7.52 ± 0.11 


-0.02 


4.542 


3.741 


2.939 


2.376 


2.020 


1.948 


1.858 


1.871 


W Sgr 


0.881 


0.109 ±0.007 


2.28 ±0.20 


8.21 ± 0.19 


-0.06 


5.412 


4.665 


3.846 


3.252 


2.952 


2.828 


2.721 


2.719 


XSgr 


0.846 


0.227 ± 0.013 


3.00 ±0.18 


7.61 ± 0.13 


-0.03 


5.299 


4.550 


3.649 


2.916 


2.626 


2.636 


2.423 


2.409 


YSgr 


0.761 


0.203 ±0.007 


2.13 ±0.29 


8.36 ± 0.30 


-0.15 


6.595 


5.738 


4.789 


4.086 


3.690 


3.609 


3.486 


3.483 


<5 Cep 


0.730 


0.073 ±0.007 


3.66 ±0.15 


7.18 ± 0.09 


-0.01 


4.602 


3.955 


3.168 


2.689 


2.378 


2.310 


2.221 


2.217 


FF Aql 


0.650 


0.204 ± 0.008 


2.81 ± 0.18 


7.76 ± 0.14 


-0.03 


6.135 


5.372 


4.488 


3.864 


3.565 


3.478 


3.378 


3.353 


T Vul 


0.647 


0.073 ± 0.008 


1.90 ±0.23 


8.61 ± 0.26 


-0.12 


6.385 


5.747 


5.015 


4.546 


4.256 


4.192 


4.114 


4.111 


RT Aur 


0.571 


0.058 ± 0.011 


2.40 ± 0.19 


8.10 ±0.17 


-0.05 


6.050 


5.468 


4.730 


4.251 


3.982 


3.924 


3.853 


3.849 
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Table 11. Multi- wavelength Period-Luminosity Relations for the LMC and MW. 









LMC 






MW 








Filter 


A [^tm] 


N 


a 


b 


<7 


N 


b 


a 


t' 


A x 


B 


0.44 


384 


-2.47 ±0.11 


15.60 ± 0.03 


0.42 


10 


-3.37 ± 0.07 


0.24 


18.97 ± 0.08 


0.51 


V 


0.55 


388 


-2.77 ± 0.08 


14.81 ± 0.02 


0.31 


10 


-4.08 ±0.04 


0.15 


18.88 ± 0.05 


0.39 


lc 


0.80 


392 


-2.97 ± 0.05 


13.93 ±0.01 


0.20 


10 


-4.77 ± 0.03 


0.12 


18.70 ± 0.04 


0.23 


.1 


1.24 


82 


-3.15 ±0.07 


13.26 ± 0.02 


0.15 


10 


-5.33 ±0.03 


0.09 


18.59 ± 0.04 


0.11 


11 


1.66 


82 


-3.23 ± 0.06 


12.90 ± 0.02 


0.12 


9 


-5.65 ±0.03 


0.09 


18.54 ± 0.03 


0.07 


K s 


2.16 


82 


-3.27 ± 0.05 


12.80 ± 0.02 


0.11 


10 


-5.70 ±0.03 


0.10 


18.50 ± 0.03 


0.04 


[3.6] 


3.55 


80 


-3.31 ± 0.05 


12.70 ± 0.02 


0.11 


10 


-5.79 ±0.03 


0.10 


18.50 ± 0.03 


0.02 


[4.5] 


4.49 


80 


-3.21 ± 0.06 


12.69 ± 0.02 


0.11 


10 


-5.77 ± 0.03 


0.10 


18.46 ± 0.03 


0.01 



Note. — The form of the PL relation is: M\ = a(logP — 1) + b. The LMC slope is used to constrain the MW 
zero-point. The apparent distance moduli (fi) are found by differencing the LMC and MW zero-points and the 
extinction (Ax) is the solution from Figure [8] 
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Table 12. Constants used in the photometric reduction. 



Channel Criterion 


u 


V 


fcorr 


zp 


Ch 1. non-saturated 


1.000 


0.080 


1.021 


280.9 ±4.1 


Ch 1. saturated 


1.000 


0.000 


1.021 


280.9 ±4.1 


Ch 2. non-saturated 


1.000 


0.025 


1.012 


179.7 ± 2.6 


Ch 2. saturated 


0.965 


-0.050 


1.012 


179.7 ± 2.6 



Note. — A source is considered saturated if 1 or more pixels 
were masked during the PRF fitting procedure; see text. 
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S Vul p = 68.464 GY Sge P = 51.36 SV Vul P = 45.0121 




PHASE PHASE PHASE 



U Car P = 38.80815 (Car P = 35.55048 T Mon P = 37.0317 




0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 

PHASE PHASE PHASE 



WZ Sgr P = 31.85 RU Set P = 19.7044 SZ Aql P = 17.14 




Fig. 1. — IRAC [3.6] and [4.5] light-curves for 37 Galactic Cepheids. The error bars corre- 
spond to the rand om photometric error and the solid line is the GLOESS interpolated curve. 
Data points from iMarengo et ajj (120101 ) are over-plotted as open triangles when available. 




Fig. I.- 



continued. 
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Fig. I. — continued. 
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Fig. 2. — Leavitt PL Relations for the HST Calibrators in the Galaxy at 3.6 (upper plot) 
and 4.5 /zm (lower plot). The data have been corrected for Lutz-Kelker-Hanson bias. The 
relations are shown as solid lines which were determined usi ng a fixed slope found from the 
LMC data and a zero-point found from the HST parallaxes ([Benedict et al.l 120071 ); the ±2<r 
boundaries are shown as dotted lines. The highly correlated magnitude residuals are plotted 
in the inset, showing that the peak-to-peak width of the IS, as defined by the HST Galactic 
Calibrators is less than 0.4 mag. 
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Fig. 3. — Leavitt PL Relations for the HST Calibrators in the Galaxy at 3.6 (upper plot) 
and 4.5 jum (lower plot). The data have not been corrected for Lutz-Kelker-Hanson bias. 
The relations are shown as solid lines which were determined usin g a fixed slope found from 
the LMC data and a zero-point found from the HST parallaxes ([Benedict et al.ll2007| ); the 
±2<T boundaries are shown as dotted lines. The highly correlated magnitude residuals are 
plotted in the inset, showing that the peak-to-peak width of the IS, as defined by the HST 
Galactic Calibrators is less than 0.4 mag. 
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Fig. 4. — Leavitt PL Relations for the Cluster Cepheids in the Galaxy at 3.6 (upper plot) 
and 4.5 jum (lower plot). The relations are shown as solid lines which were determined using 
a fixed slope found from the LMC data. The ±2er boundaries are shown as dotted lines. The 
highly correlated magnitude residuals are plotted in the inset, showing that the peak-to-peak 
width of the IS, as defined by the Cluster Galactic Calibrators is less than 0.8 mag. 
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Fig. 5. — Leavitt PL Relations for the IRSB Cepheids in the Galaxy at 3.6 (upper plot) and 
4.5 fim (lower plot). The relations are shown as solid lines which were determined using a 
fixed slope found from the LMC data. The ±2cr boundaries are shown as dotted lines. The 
highly correlated magnitude residuals are plotted in the inset, showing that the peak-to-peak 
width of the IS, as defined by the IRSB Galactic Calibrators is less than 0.6 mag. 
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Fig. 6. — The Galactic Period-Color Relation. The solid circles are the de-reddened colors 
for the 37 Cepheids in this study; see Tables Q]& HI The so lid line i s the fit to all the data 
except Y Oph. The broken lines are theoretical trends from iNgeowl (120121 ). 



o° ° 

i i i i L 



1 1.5 
Log Period [days] 



Fig. 7. — The Period Luminosity laws for the B, V, Ic, J, H, K s , [3.6] and [4.5] bands found 
for LMC data (open gray circles) and the Galaxy (filled black dots). Solid lines are the PL 
relations for each band corrected for LMC distance and extinction. 
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Fig. 8. — Standard extinction curve fit to multi-wavelength apparent distance moduli to the 
LMC. The fitted extinction model yields a true distance modulus (intercept) of 18.48 ±0.03 
and an average color-excess (slope) of E(B — V) = 0.12 ±0.01 mag. The K s and [4.5] bands 
were not included in the fit due to the effect of CO in those bands. 
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Fig. 9. — De-reddened color versus Period. The red extent of each Cepheid is well defined 
near 0.01 mag. The blue extent increases with period due to CO molecular absorption in 
Cepheids that reach intrinsically cooler temperatures. 
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Fig. 10. — Left: Synthetic spectra showing the onset of CO molecular absorption in the 
IRAC Channel 2 bandpass. Right: The synthetic color, normalized to zero at 7000K, for 
the generated model spectra. The color trends to the blue at cooler temperatures, a result 
of the CO molecular absorption. 
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Fig. 11. — The pixel phase correction (PPC) as a function of pixel phase for a sample of 
Galactic Cepheids. The black points are for non-saturated sources and the gray points are 
for saturated sources; see text. Top left: The Channel 1 data as output from the MOPEX 
PRF fitting algorithm. Top right: The Channel 1 data after the pixel phase correction has 
been applied. Bottom panels: same as top, but for Channel 2 data. 
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Fig. 12. — Dither time-series of PSF subtracted images for one epoch of RT Aur showing the 
9 dither positions (Channel 1). RT Aur has been fit using MOPEX and the model fit was 
subtracted leaving the residual. In these frames (ordered 1-9) it is possible to see the latent 
image of the star trail as the telescope slewed into position at time-stamp 1. The trail is 
nearly dissipated by time-stamp 9; but new latent images are formed at the location of each 
the previous dithers. By time-stamp 9 the previous 8 latent images reveal the full Reuleaux 
dither pattern as a faint after glowing latent triangle. 
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Fig. 13. — Image persistence mitigation for Channel 1. Top left: Original BCD image of 
the Galactic Cepheid RT Aur. Top middle: Residual image after preliminary PRF fitting. 
Top right: Persistence map created from weighted average of nine dithered residual images. 
Bottom left: Original BCD with background/persistence map subtracted. Bottom middle: 
The uncertainty map with background pixels masked and only the core region used for PRF 
fitting. Bottom right: Final residual image after persistence map was subtracted. The same 
procedure was followed for Channel 2 although the latency effect was negligible. 



-47- 



REFERENCES 

Barnes, III, T. G., Fernley, J. A., Frueh, M. L., Navas, J. G., Moffett, T. J., & Skillen, I. 
1997, PASP, 109, 645 

Benedict, G. F., et al. 2007, AJ, 133, 1810 

Berdnikov, L. N. 2008, VizieR Online Data Catalog, 2285, 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Castelli, F., & Kurucz, R. L. 2003, in IAU Symposium, Vol. 210, Modelling of Stellar Atmo- 
spheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, 20P 

Chapman, N. L., Mundy, L. G., Lai, S.-P, & Evans, II, N. J. 2009, ApJ, 690, 496 

Contreras, C, et al. 2010, AJ, 139, 519 

Fazio, G. G., et al. 2004, ApJS, 154, 10 

Fernie, J. D., Evans, N. R., Beattie, B., & Seager, S. 1995, Information Bulletin on Variable 
Stars, 4148, 1 

Flaherty, K. M., Pipher, J. L., Megeath, S. T., Winston, E. M., Gutermuth, R. A., Muzerolle, 
J., Allen, L. E., & Fazio, G. G. 2007, ApJ, 663, 1069 

Folatelli, G. 2009, in Revista Mexicana de Astronomia y Astrofisica, vol. 27, Vol. 35, Revista 
Mexicana de Astronomia y Astrofisica Conference Series, 219-220 

Fouque, P., et al. 2007, A&A, 476, 73 

Freedman, W., Madore, B., Mager, V., Persson, E., Rigby, J., & Sturch, L. 2008, in Spitzer 
Proposal ID #60010, 60010 

Freedman, W. L., & Madore, B. F. 2010, ARA&A, 48, 673 

Freedman, W. L., et al. 2011, AJ, 142, 192 

-. 2012, ApJ, Accepted 

Hanson, R. B. 1979, MNRAS, 186, 875 

Indebetouw, R., et al. 2005, ApJ, 619, 931 

Kovtyukh, V. V., Soubiran, C, Luck, R. E., Turner, D. G., Belik, S. I., Andrievsky, S. M., 
& Chekhonadskikh, F. A. 2008, MNRAS, 389, 1336 



-48- 

Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. Kurucz CD-ROM 
No. 13. Cambridge, Mass.: Smithsonian Astrophysical Observatory, 1993., 13 

Laney, C. D., Joner, M. D., & Pietrzyhski, G. 2012, MNRAS, 419, 1637 

Laney, C. D., & Stobie, R. S. 1992, A&AS, 93, 93 

Loredo, T. J. 2007, in Astronomical Society of the Pacific Conference Series, Vol. 371, Sta- 
tistical Challenges in Modern Astronomy IV, ed. G. J. Babu & E. D. Feigelson, 121 

Lutz, T. E., & Kelker, D. H. 1973, PASP, 85, 573 

Majaess, D., Turner, D., & Gieren, W. 2012a, ApJ, 747, 145 

Majaess, D., Turner, D., Gieren, W., Balam, D., & Lane, D. 2012b, ApJ, 748, L9 

Majaess, D., et al. 2011, ApJ, 741, L27 

Makovoz, D., & Marleau, F. R. 2005, PASP, 117, 1113 

Marengo, M., Evans, N. R., Barmby, P., Bono, G., Welch, D. L., & Romaniello, M. 2010, 
ApJ, 709, 120 

Matthews, J. M., Gieren, W. P., Mermilliod, J.-C, & Welch, D. L. 1995, AJ, 110, 2280 
Monson, A. J., & Pierce, M. J. 2011, ApJS, 193, 12 
Ngeow, C.-C. 2012, ApJ, 747, 50 

Ngeow, C.-C, Marconi, M., Musella, I., Cignoni, M., & Kanbur, S. M. 2012, ApJ, 745, 104 

Nishiyama, S., Tamura, M., Hatano, H., Kato, D., Tanabe, T., Sugitani, K., & Nagata, T. 
2009, ApJ, 696, 1407 

Persson, S. E., Madore, B. F., Krzemihski, W., Freedman, W. L., Roth, M., & Murphy, 
D. C. 2004, AJ, 128, 2239 

Reach, W. T., et al. 2005, PASP, 117, 978 

Samus, N. N., Durlevich, O. V., & et al. 2009, VizieR Online Data Catalog, 1, 2025 

Sbordone, L. 2005, Memorie della Societa Astronomica Italiana Supplementi, 8, 61 

Sbordone, L., Bonifacio, P., Castelli, F., & Kurucz, R. L. 2004, Memorie della Societa 
Astronomica Italiana Supplementi, 5, 93 



-49- 



Scowcroft, V., Freedman, W. L., Madore, B. F., Monson, A. J., Persson, S. E., Seibert, M., 
Rigby, J. R., & Sturch, L. 2011, ApJ, 743, 76 & erratum 

Scowcroft, V., Freedman, W. L., Madore, B. F., Monson, A. J., Rigby, J., Persson, S. E., & 
Sturch, L. 2012, In Prep 

Smith, H. 2003, MNRAS, 338, 891 

Storm, J., Gieren, W., Fouque, P., Barnes, T. G., Soszyhski, I., Pietrzyhski, G., Nardetto, 
N., & Queloz, D. 2011a, A&A, 534, A95 

Storm, J., et al. 2011b, A&A, 534, A94 

Tammann, G. A., Sandage, A., & Reindl, B. 2003, A&A, 404, 423 

Turner, D. G. 2010, Ap&SS, 326, 219 

— . 2011, Rev. Mexicana Astron. Astrofis., 47, 127 

Turner, D. G., & Burke, J. F. 2002, AJ, 124, 2931 

Udalski, A., Soszynski, I., Szymanski, M., Kubiak, M., Pietrzynski, G., Wozniak, P., & 
Zebrun, K. 1999, Acta Astron., 49, 223 

Walker, A. R. 2011, Ap&SS, 746 

Welch, D. L., Wieland, F., McAlary, C. W., McGonegal, R., Madore, B. F., McLaren, R. A., 
& Neugebauer, G. 1984, ApJS, 54, 547 

Windmark, F., Lindegren, L., & Hobbs, D. 2011, A&A, 530, A76 



This preprint was prepared with the A AS IATgX macros v5.2. 



